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ABSTRACT 


The need for a proper theoretical explanation of the residual 
film, reported Ly various workers, has led to the development of 
microcontinuum theory of lubrication. Even after a decade this 
theory seems to be in a preliminary state. Scores of lubrication 
problems are yet to be analysed from this view point. For example, 
little work has appeared on dynamically loaded bearings and porous 
bearings from the microoont inuum view point. Moreover, no attempt 
has so far been made to study the impact of the micropolarity 
on finite dimensional bearings and bearings with rough surfaces. 

Thus, in this thesis, an endeavour is made to study these 
lubrication problems using the micropolar fluid theory. The thesis 
is divided into IX Chapters. 

Chapter I is of introductory type, wherein are discussed 
important experimental evidences, which establish the existence of 
the residual film (rheological abnormalities) . It also emphasizes 
the fact that the origin of the so called rheological abnormalities 
may lie in the micro-mcchanical motion of the fluid element. Thus 
this Chapter lays the foundation for the present work. 

In Chapter II, the basic equations for the micropolar fluid 
are given. Attention is restricted to a three dimensional laminar 
flow of an incompressible micropolar fluid in a finite slider type 
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"bearing, Application of usual lubrication assumptions and a subsequent 
detailed order of magnitude analysis reduces the governing equations 
to two systems of coupled ordinary differential equations to be solved 
for velocity distributions and microrotation velocities (no paper, 
published so fab, gives such a detailed order of magnitude analysis 
for the three dimensional problem) , Generalized Reynolds equation 
is then obtained in a conventional manner by integrating the continuity 
equation. Also given in this Chapter are expressions for flow fluxes 
and shear stresses. 

In Chapter III, the generalised Reynolds equation, obtained 
in Chapter II, is applied to analyse the problem of two dimensional 
squeeze film of a ball in a spherical seat and to some three dimensional 
r.on-cyclic squeeze films, in an effort to study the effects of rigid 
particle additives for the three dimensional raicropolarity model. 

In Chapter IV, dynamically loaded short journal bearings are 
studied, Reynolds equation for the general cas^, \7hero both wedge 
and Squeeze films are active, is derived. Detailed consideration is 
given to the dynamic behaviour of squeeze film in a short journal 
bearing subject to a sinusoidal loading with no journal rotation. 

Chapter V, also gives an analysis of dynamically loaded hearings* 
The problem considered is that of the roller bearing in combined rolling, 
sliding and normal motion, under cavitation boundary conditions. 



In Chapter 71, the micropolar theory is applied, to the analysis 
of hearings with rough surfaces. Generalized forms of Reynolds 
equation are derived for three cases of roughness, namely, longitudinal, 
transverse and isotropic, using the stochastic approach. These 
equations are subsequently applied to the problem of infinitely 
long journal bearing whose stationary surface (bearing) is assumed 
to be rough, in an order to study the interaction of micropolarity 
with surface asperities, using the half Sommerfeld boundary conditions. 

In Chapter YII, the generalized Reynolds equations, obtained 
in Chapter 71, are applied to study the squeezing between parallel 
plates. Three geometries, namely, infinitely long parallel plates, 
circular plates and finite rectangular plates are considered. 

In Chapter Till, the problem of synovial joints lubrication is 
studied, by approximating the joints by a spherical bearing, the 
cartilage surface as porous (but rigid) and the lubricating synovial 
fluid as micropolar. A modified form of Reynolds equation is derived 
using the slip boundary conditions. 

Chapter IX is of conclusive nature wherein are presented the 
main conclusions of the work done in the thesis and the conclusions 
are compared with the various experimental evidences. The one overall 
conclusion that emerges is the existence of increased effective 

v 

viscosity in the proximity of a solid surface. This increase is further 
enhanced when the surfaces are rough. 



vii 


A part of the work presented in this thesis has been accepted 
for publication in the form of following research papers. 

(1) Theoretical Effects of Rigid Particle Additives in Ron-Cyclic 
Squeeze Films, *TRAKS .ASMS, J .Lub.Tech’ . 

(2) Dynamic Loading in Micropolar Fluid Lubricated Short Bearings. 
1 J .Mech .Engng .Sci . , Instn . of Mech . Engrs' . 



viii 


LIST OP PAPERS 

During the research tenure (December 1977 to March 198 l) 
thirty ( 30 ) research problems were finalized. These are listed in 
two categories. 

(1) Category A - Related to the thesis 

(2) Category B - Not related to the thesis 

The details of each category is as follows 

CATEGORY A - WORK RELATED TO THE THESIS 

The work presented in the thesis has been accept ed/communicated 
for publication in the forms of the following research papers t 

(i) Papers Accepted 

(1) Theoretical Effects of Rigid Particle Additives in Ron-Cyclic 
Squeeze Films . 

TRANS .ASME, J .Lub .T ech . , ( IN PRESS) . 

(2) Dynamic Loading in Micro polar Fluid Lubricated Short Bearings. 
J.Mech.Engng.Sci.,Instn.Mech.Engrs. (IN PRESS). 

(ii) Papers Communicated 

( 3 ) Three Dimensional Reynolds Equation for Micropolar Fluid 
Lubricated Bearings. 

( 4 ) Lubrication of Roller Bearing in Combined Rolling, Sliding 
and Normal Motion With Additives • 

( 5 ) Lubrication of Rough Surfaces - A Microcontinuum Analysis. 

( 6 ) Microcontinuum Analysis of Squeeze Films Between Rough Surfaces 

( 7 ) Micropolar Squeeze Films Between Rough Rectangular Plates. 

CO 

( 8 ) Lubrication of Human Joints - A Microfat inuum Approach. 



CATEGORY B-WOEK ROT RELATED TO THE THESIS 


ix 


Apart from the work presented, in the thesis a number 
of other research problems were finalized.. Details are as follows r 


(i) P apers Published. 

(9) The Effect of Additives in the Lubricant of a Composite 
Bearing with an Inclined Stepped Surface. 

W e ar,Yol. 66 ,No . 1 , 1981, p.17-26 . 

( 10 ) Effect of Viscosity Variation Due to Lubricant Additives 
in J ournal Bearings . 

W e ar,Yol. 66 ,No .1, 1981 . 

(11) Couple Stresses in Journal Bearing Lubricants and the Effect 
of Cavitation. 

Wear, Yol. 67, Ho. 1,1981, p.15-24. 

(1 2 ) Couple Stresses in the Lubrication of Rolling Contact 
Bearings Considering Cavitation. 

Wear,Yol.67,Eo.l, 1981, p.85-97 « 

(ii) Papers Accepted 

( 13 ) Lubrication of a Cylinder on a Plane With a Non-Newtonian 
Fluid Considering Cavitation. 

TRANS .ASME,J.Lub. Tech. (IN PRESS). 

( 14 ) Non-Newtonian Squeeze Filins in Journal Bearings. 

Wear. ( IN PRESS) . 

( 15 ) Couple Stresses in the Elastohydrodynamic Film in Roller 
Bearings . 

Wear. (IN PRESS). 

(16) Non-Newtonian Squeeze Films in Spherical Bearings. 

Wear . (IN PRESS) . 

(iii) Papers Comiomicated 

(17) Non-Newtonian Power-Law Fluid Lubrication of Lightly Loaded 
Cylinders With Normal and Rolling Motion. 



X 


(18) Cyclic Squeeze Film in Porous Bearing Using Velocity Slip 
Conditions. 

( 19 ) A Microcontinuum Analysis of the Self Propulsion of the 
Spermatozoa in the Cervical Canal. 

(20) Elastohydrodynamic Lubrication of Circular Plate Thrust 
Bearing With Power Law Lubricants. 

(21) Viscosity Variation in Journal Bearing Lubricant With Additives 
Considering Cavitation. 

(22) Effect of Additives on Finite Step Bearing Lubricant. 

(23) Non-Newtonian Lubrication Theory for Rough Surfaces ; 
Application to Rigid and Elastic Rollers. 

(24) Effects of Surface Roughness and Additives in Lubrication : 
Generalized Reynolds Equation and Its Application to 
Elastohydrodynamic Film. 

(25) Non-Newtonian Theory for Rough Squeeze Films. 

(26) Dynamically Loaded Rough Journal Bearings . 

(27) Radial and Frictional Forces in Misaligned Radial Face Seals 
With Non-Newtonian Fluids . 

(28) Hydrostatic Pressure Effects in Misaligned Radial Face 
Seals With Non-Newtonian Fluids. 

(29) Elastohydrostatic Forces in Misaligned P-adial Face Seals. 

(30) Micropolar Squeeze Films in Porous Hemispherical Bearings. 



CONTENTS 


xi 


CHAPTER 

CERTIFICATE 

ACKNOWLEDGEMENTS 

ABSTRACT 

LIST OF PAPERS 

NOMENCLATURE 


I s GENERAL INTRODUCTION 1 



1.1 

Why lubrication 

1 


1.2 

Early days in lubrication 

2 


1.3 

Classification of lubrication 

4 


1.4 

Lubrication and rheology 

5 


1.5 

Experimental evidences 

7 


1.6 

Developments in theoretical rheology 

13 


1.7 

Microcontinuum analysis of fluid flows 

17 

II * 

BASIC EQUATIONS AND GENERALIZED REYNOLDS EQUATION FOR 
IICROPO'LAR FLUIDS 

21 


2.1 

Introduct ion 

21 


2.2 

Balance equations for micropolar fluids 

22 


2.3 

Lubrication problem-As sumptions and basic equations 25 

Ills 

TWO 

AND THREE DIMENSIONAL NON-CYCLIC SQUEEZE FILMS 

37 


3.1 

Introduction 

37 


3.2 

Spherical bearings 

38 


3.3 

Sliders and rectangular plates 

42 


3.4 

Elliptical plate® 

44 


3.5 

Miscellaneous configurations 

45 


3.6 

Results and discussion 

50 

17 s 

DYNAMICALLY LOADED SHORT JOURNAL BEARINGS 

65 


4.1 

Introduction 

:.&5 


4.2 

The problem and analysis 

67 


4-3 

Results and discussion 

75 

7 i 

LUBRICATION OF ROLLER BEARINGS IN COMBINED ROLLING, 
SLIDING AND NORMAL MOTION 

86 


5.1 

Introduct ion 

86 


5.2 

The problem and theoretical analysis 

87 


5.3 

Results and discussion 

96 


i 

ii 
iv 

viii 

xiii 



xii 


VI % LUBRICATION THEORY FOR ROUGH SURFACES AND ITS 

APPLICATION TO A JOURNAL BEARING 108 

6 .1 Introduction 108 

6 .2 Modified. Reynolds equation for rough surfaces 110 

6.3 Bearing characteristics 116 

6 .4 Modified Reynolds equation applied to a journal 

hearing 118 

6.5 Non-dimensional forms 123 

6 .6 Results and discussion 126 

VII : SQUEEZE FILMS BETWEEN ROUGH SURFACES 138 

7 .1 Introduction 138 

7.2 Mathematical analysis 141 

7.3 Squeeze films Between rough infinitely long 

parallel plat es 143 

7.4 Squeeze film Between rough circular plates 147 

7.5 Squeezing Between rough rectangular plates 151 

7.6 Numerical results 157 

7.7 Results and discussion 158 

VIII s LUBRICATION OF HUMAN JOINTS 170 

8 .1 Introduction 170 

8.2 Anatomy of a synovia 1 joint 176 

8.3 Theoretical formulation-As sumptions and equations Yf6 

8 .4 Normal motion of human joint 186 

8.5 Results and discussion 191 

IX j SUMMARY AND CONCLUSIONS 207 

REFERENCES 217 



xiii 


a 

a 


1 

2 


A 

A 

b 

b 

B 

B 


1 


2 


1 


c 


J R 


D 


E(f) 

f(h g ) 

£(N,a f h) 


F 

F 


1,2 


F R, ’ F R r 


R 


HOMMCLATURE 

eccentricity for spherical hearing 

x-line ns ion of rectangular plate (or slider hearing) 

area of rectangular plate 

major axis of elliptical plate 

half total range of random film thickness variable 

z-dimension of rectangular plate (or slider hearing) 

roughness parameter 

minor axis of elliptical plate 

quantity, defined by eqn. (5*23) 

hearing clearance 

quantity defined by eqn. (5*24) 

coefficient of friction parameter 

diameter of the journal hearing 

eccentricity of the journal hearing 

expected or mean value of f 

roughness distribution function 

defined by eqn. (2.52) fat Chapters' III, IV and V, and defined 
by eqn. (2,55) for Chapters Vi and VII 

frictional force 

frictional forces on the plane and cylindrical surface, 
respectively 

frictional force ratio parameters 
frictional drag ratio parameter 


t 



xiv 


p* 

P(N # L,H) 

]? 1 (N,L f B) 

P 2 (N,L,H n ,B) 

f 3 (k,l,b) 

g(N,A,h) 


g(n,l,h) 

g(n,l,h,b) 

h 


h 


hn 


h 


n 


' m 


h 


h. 


'n. 


H tf u ti U 

7 n 7 nr s 7 1 
3 


J 

£ 


£ 


1 


L 



L 


o 


body force per unit mass 
defined by eqn. (3.14) 
defined by eqn. (7*15) 
defined by eqn. (7 .19) 
defined by eq.ru (7 .27) 
defined by eqn. (7.28) 

defined by eqn. (2*6l) for Chapters III, TV and V ? 
and by eqn. (2.64) for Chapters VI and VII 

defined by eqn. (5*20) 

defined by eqn. (6.55) 

film thickness 

minimum/f inal film thickness 
initial film thickness 
nominal (average) film thickness 
film thickness corresponding to x = x ffi 
random part of the film thickness 
initial (nominal) film thickness 

non-dimensional form of h,h^,h m ,h s ,h^ | respectively 
microinertia constant 
defined by eqn. (5*20) 

characteristic material length defined by eqn. (2.18) 
half x-dimension of infinitely long parallel plate 
characteristic length ratio parameter 
body couple 

length of the short journal bearing (in axial direction) 





XV* 

m 

F 

£ 



couple stress tensor 


N 

coupling number, defined by eqn. (2 .18) 


n 

revolutions per unit time 


P 

hydrodynamic pressure 


P 

non-dimensional pressure 


P o 

unit loading, — — 

0 


q 

velocity ratio parameter, defined by eqn. (5*6) 


q. 

flow flux (in vector form) 


ws 

fluxes along x,z and 0 directions, respectively 


r 

radial co-ordinate 


r l ’ r 2 

radii of circular plates 


R 

radius of journal/cylinder/sphere 


R 1 

larger radius of conical seat 


r 2 

smaller radius of conical seat 


Re 

modified Reynolds number 


Rte 

defined by eqn. ( 2 . 25 ) 


S 

Sommerfeld number 


t 

time 


^k£ 

stress tensor 


T 

non-dimensional time 


% 

time ratio parameter 


u,v,w 

velocity components 


U 

velocity of the journal 




xv i 


UpUg 

T 

Y 
-»■ 
t b 

V 

p 

II 

W(t) 

w. 


w .w 

o* 7r 


i 


s 

x,y,z 


x,x 


m 


P 


0 
<fr 

Ip 

a , 6 ,Y,X 

“l 


velocity of the plane and cylindrical surface, respectively 
normal relative velocity of the cylindrical surfaces 
velocity vector 

velocity vector on the solid boundary 
velocity of the fluid towards porous region 
load capacity- 

load capacity as a function of time 
amplitude of the applied load 

load components along and perpendicular to the line 
of centers, respectively 

load capacity ratio parameter 

cartesian co-ordinates 

non-dimensional form of x and x m , respectively 
angular co-ordinate (latitude angle) 
angle at which p = o 

angle between load line and line of centers 
porosity of the cartilage surface 
non-dimensional porosity parameter 
non-dimensional slip parameter 
angle between the vertical line and load line 
viscosity coefficients for a micropolar fluid 
semi-vertical angle of cone 

angle subtended by the radii r^ and , at the center 
of circular sector 


e 


eccentricity ratio 



C4* 


xvii 


>* e i 


£ 

k £ r 
X , y 


B 


V v 2 ,v 3 


c 

0 ) 


CD 

P 

P 

a 

T 


T 


Z 


n 


e 

n 

c 

A 


eccentricity ratios corresponding to the final and 
initial film thickness 

alternating tensor 

Newtonian viscosity coefficients 

dynamic viscosity coefficient 

coefficient of friction 

microrotation vector 

micro rot at ion vector on the "boundary 

components of the microrotation vector 

aspect ratio 

critical aspect ratio 

angular velocity of the journal 

d^ 

frequency of the rotation of the applied load 5 

frequency of oscillation 
mass density 
standard deviation 
shear stress 

shear stresses along x-direotion and z-direction, respectively 
angular velocity of the cylinder 
random variable 
thermodynamic pressure 

longitudinal angle (used for spherical polar co-ordinates) 
slip parameter 

•A bar above a variable denotes the corresponding non-dimensional 


variable 



CHAPTER I 


GENERAL INTRODUCTION 


1.! WHY LUBRICATION 

Lubrication plays a most vital role in our great and complex 
civilization. To estimate the importance of this role one need only 
consider that every moving part of every machine is subjected to 
friction and wear. Throughout the centuries from ox-cart axles to 
grinding wheel spindles, one of man’s most persistent problems has 
centered around reduction and control of friction and wear. Friction 
consumes and wastes energy. Wear causes changes in dimensions and 
eventual breakdown of the machine element and the entire machine and 
all that depends upon it. 

Today the constituent parts of modern high speed engines are 
frequently subjected to a very high rate of we ar, which can lead to 
serious operational disturbances, if not to the total deterioration 
and destruction of engines* The need for continuous engine operation 
today necessitates the timely removal and replacement of worn engine 
parts, which can be a costly operation. A nation’s economy can be 
seriously impaired through excessive losses of engine and machine 
parts due to wear, 'so that an intensive campaign against friction and 
wear is in the best of interest of increased production and national 
efficiency. 

So whether one is concerned with reducing the waste of world’s 
production of energy, conserving a nation’s critical natural resources. 



maintaining production schedules in a plant, creating a larger 
margin of profit in plant operation, or even making the family 
automobile run longer and better, one must be interested in the 
study of lubrication, 

1.2 JEARLf DAYS IN LUBRICATION 

Lubrication, in one sense, is as old in civilization as the 
wheel and axles. When a tomb was opened in Egypt, some years ago, 
one of the chariots still had some of the original lubricant on the 
axle. This was analyzed and found to be sticky and slightly 
greasy [ 1 ] • It contained road dirt such as quartz sand, compounds 
of aluminum, iron, and lime. It had a melting point of ^O 0 ? 1 , which 
suggests that it might have been mutton or beef tallow, either of which 
would have proved suitable for axle lubrication in that warm country. 

Until comparatively recent years all lubricants were largely 
of animal, vegetable, or marine origin : mutton tallow, lard, gooze 
grease, fish oils, castor oil, cottonseed oil and other vegetable 
oils, etc. The use of mineral oils dates back principally a little < 

i 

over a hundred years . > 

f 

Pioneering work in the field of lubrication dates back to the ! 
year 1847* when Adams [ 2 ] developed and patented several rather good j 
designs for railway axle bearings* Many of the plans show a remarkable J 
comprehension of the fundamentals of proper bearing design. 

The new science of lubrication, however, based on hydrodynamics : 
and oil flow, was not recognized as such until 1883* when Petrov [ 3 i 
made the first significant attempt to analyze theoretically the friction 



effect of film lubrication. It was then recognized that the film rather 
than the hearing material could he a prime consideration. 

In the same year that Petrov published his theory, unexpected 
experimental results were reported hy Tower [ 4 ] . He demonstrated 
that a loaded, oil-lubricated journal hearing is subject to local 
pressures substantially higher than its mean pressure. Three years 
later Reynolds [ 5 ] , unaware of Petrov’s theory, was able to explain 
Tower’s results. Prom the properties of a thin laminar fluid film 
he derived and published not only the descriptive differential equation, 
which even today bears his name, but also certain solutions to this 
equation that agree well with the experimental measurements of 
Tower [4 ] . Several years later in 1904 , Sommerfeld [6] , provided 
elegant theoretical extensions to the journal bearing lubrication 
problem • 

Of the thousands of papers on film lubrication which have 
appeared since 1886, in the wake of Reynolds classical work on lubrication 
theory, few have offered modification of the theory. In general, papers 
on film lubrication have extended the application of the fundamentals, 
as set forth by Reynolds, to a variety of bearing shapes with steady 
and transient films of variable viscosity, density and temperature. 

In 1928, Boswall [ 7 ] summarized the theory of film lubrioation to 
that date. By now, however, much additional information about lubricating 
films has been developed. A proper historical review of 1 *be same 
would be a Herculean task and beyond the scope of the present thesis. 
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Inspite of the rapid growth of the lubrication theory in the 
"fcii 

20 century, the picture is far from complete. A perfect correlation 
between theory and practice has not yet been obtained. This is fundamental! 
because the basic assumptions frequently circumscribe a rather tight circle 
within which the applicability of the theory is valid . 

1.3 CLASSIFICATION OF LUBRICATION 

In surveying the entire field of lubrication with its many 
vital physical and chemical manifestations, the one phenomenon that 
occurs most frequently and persistently and which probably has by far the 
most value is the creation of a load-carrying fluid film. For the builders 
and users of machinery and mechanism of various kinds, development of 
this fluid film is quite a natural action and serves them obediently 
and tirelessly if given even the slightest encouragement. It is obvious 
that much of our industrial life would grind to a screeching halt without 
this lubrication phenomenon. 

The thickness of the fluid film not only determines between 
the various kinds of lubrication but also the regime of applicability 
of the theory. At one extreme is the thick film lubrication, where 

the bearing surfaces are completely separated by a continuous fluid 

-3 

film of the order of 10 cm. and' the bearing characteristics mainly 
depend on the viscosity of the lubricant in bulk. At the other extreme 
is the boundary lubrication, where a fluid film only a few molecular 
layer thick is left between the bearing surfaoes • Between these two 
extremes, lies, the region of thin film or mixed lubrication. The fluid 
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film is very thin, of the order of 10"^ cm. so that the load is partly 
supported "by the roughness asperities and partly "by the fluid film. The 
region of thin film lubrication is an indeterminate state in which the 
strict laws of hydrodynamic lubrication are no longer valid. Some of the 
most exciting research is being done in this field. It is actually 
a meeting place for a chemist and an engineer as well as a metallurgist 
and a mathematician. It is in this region that the rheology plays an 
important role. It is in this region that the proximity of a solid 
surface drastically modifies the properties of the lubricant. It is 
also an area of great endeavour, where the experimental facts are to be 
correlated with theory. 

1.4. LUBRICATION AID RHEOLOGY 

Interaction between the fields of lubrication and rheology has 
become increasingly important within the past few years with the growing 
awareness of the complexity of the lubricant reaotion to severe and ; 

transient stress conditions, the continued demand for well lubricated, 
high speed, highly loaded machine elements, extensive use of synthetic 
lubricants and use of additives. j 

f 

The general definition of rheology is the science that deals ! 

I" 

with the deformation and flow of materials. Similarly, one may state, ; 

that the primary objective of lubrication is to separate rubbing surfaoes ' 

t 

% 

by a layer of lubricant to prevent excessive friction and wear. It is cleg 
that the lubricant during the lubrication process will be subjected to a 
quite high shearing stresses and velocity gradients. The rheological 
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properties of the lubricants are, therefore, likely to have an important 
effect on its action. 

In early lubrication work, only the Newtonian viscosity of the 
lubricant was important. However the applied stresses have increased 
and the range of environmentation conditions, under which the lubricant 
should work, have become more extreme. The fact that many mineral 
base oils also manifest non-Newtonian behaviour under the severe shear 
conditions they often encounter in practice also support this point 
of view, 

Non-Newtonian behaviour is also almost invariably observed in 
various lubrication systems, when additives are used. These additives 
are put into the lubricant for a variety of purposes and do a great f 

deal to improve the lubricating oils which nature and the refiner produce* 
Modern lubricating-oil additives, based on years of scientific research, 
designed to meet the extreme demands of modern machines and for high 
performance ratings, under actual working conditions, have become 
indispensable in many applications. 

The ever girowing trend to transmit greater power through 
mechanisms smaller in size and weight imposes an increasing burden | 

on the lubricant. New combustion problems, greater surface loadings, 

1 

wider range of operating temperatures and often greater stalling or sliding! 

! 

I 

speeds in bearings, all subject the lubricants to higher operating j 

| 

temperatures, loads and overall abuses. ' , 

I 

These examples show just a few areas where a knowledge of the 
rheological properties of the lubricant is likely to be important. Nor 
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a proper comprehension, one mast look for the experimental evidences 

been 

where such properties have "been found, or have/sought for and found 
lacking. Hence in the next section a brief summary of such experimental 
studies is given. 

1.5 EKPEdJMEhTAL EVIDENCES 

The manner in which a liquid in contact with a solid surface 
is influenced hy the proximity of that surface has been the subject 
of intensive study and a matter of controversy. It appears obvious 
that the molecules of the liquid in the immediate vicinity or those 
adhering to the surface, or of the adsorbed layers on the surface 
must necessarily behave different from the adjacent molecules of 
the liquid in bulk. However the nature of this surface influence 
and the depth to which it penetrates is not yet properly understood 
and the experimental results have given rise to widely divergent 
opinions. It is generally accepted that the surface zone of a 
liquid is not merely a monomolecular layer but a raultimolecular 
layer [8,9 ] * i*©* the molecular orientation extends effectively 

i 

to many molecular lengths. I 

I 

From a practical point of view, however, the important question | 

j 

is not whether deep molecular orientation exists, but whether it gives : 

i 

rise to rheological abnormalities in the neighbourhood of a solid surface 
i.e. - whether a fluid in the immediate vicinity of a solid surface, devej 
properties significantly different from the adjacent molecules In bulk. 
The experimental evidences are viewed from this angle. ! 
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Kings "bury L 10 j seems to be among the earliest investigators 
of the question of existence of rheological abnormalities in proximity 
of a solid surface. He observed an enhanced viscosity in that part of the 
fluid which is in the region of attraction of the surface molecules of the 
metals. This enhancement may range from a small increase in viscosity 

_7 

to rigidity and the depth of the affected zone ranges from 10 mm to 

-3 

10 mm. 

Hardy and Nottage [ 11-13 ] offer proofs that the liquid film 
must at least be of 0.007 to 0.010 mm in thickness, if any of the liquid 
is to be beyond the range of surface influence. v/illsoR and Barnard [ 14 ] 
found that capillaries as large as 0.3 mm diameter would gradually clog 
up under the flow of oils and would close completely if a small percentage 
of oleic or stearic acid were added to the oil. This clogging was 
attributed to the gradual building up of adsorbed films on the walls of 
the capillaries. Bulkley [ 15 ] in a carefully conducted series of 
experiments with capillaries of approximately the same size as those 
used by Willson and Barnard [ 1/f.J , was however able to show that the 
clogging was due to the fact that the oil had not been thoroughly 
filtered. Prom farther experiments with highly filtered oils, Bulkley 
concluded that there is no rigidity in clean liquid at distances of 
0.3 x 10" 4 mm from a solid boundary. This, although in accordance 
with the two experiments of Bastow & Bowden [l6,17] , who detected no 
evidence of non-Newt onian behaviour upto a distance of 10"" 4 mm, is 
in c cntradiction to Griffiths [ 18 ] , who claimed to have found evidence 
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of rigidity in hydrocarbons of low molecular weight with gap between 

-3 -3 

10 mm to 2 x 10 mm. Results similar to this have also been obtained 

by Tausz and Szekely [19] , who found that a film of lubricating 

011 under a steel ball dipped into mercury ’was remarkably stable. 

Non~polar oils gave a film about 5000 A (5 x 10~ 5 mm) which broke only 
after two days, films of polar oils maintained themselves for as long 

as sixty five days. Terzaghi [ 20 J found the viscosity of water to be 

„ 5 

many times the bulk value in the passage between glass plates 10 mm 
apart. Macaulay [ 21 ] also found a similar increase in viscosity. 

Of great importance and interest is the work of Needs [ 22 ] 
in 1940, who squeezed fluids between two steel plates and found that 
a stable residual film was formed. Needs’ critics expended most of 
their energy on deciding that the fluids used were contaminated with 
solids. Needs’ defence partially allays this situation at least for 
olive and castor oils • Though one must view Needs’ work with some 
reservation, it appears very reasonable that an alteration in the 
liquid film occur ed making it substantially more viscous than in 
the bulk. 

In 1949 Henniker [ 8 ] published an extensive 
review paper which included 174 references, largly experimental ones, 
in support of his conclusion that the surface zone of a liquid is not 
merely a monomolecular layer, but a multimolecular layer. More recently, 
but much more briefly, Adamson [ 2J j in i960 admitted that although 
all the evidence of deep surface orientation is more or less circumstantial, 
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nevertheless deep molecular orientation does ocour, probably as a result 
of perturbations passed on from molecule to molecule. 

•An important class of problems where the evidence of a greatly 
enhanced viscosity in thin films close to a solid surface is strongest, 
is when small proportion of a surface active compound is added to the 
fluid. In fact, some maintenance practices recommend the addition of a 
small amount of used oil to new oil that may bo placed in a lubrication 
system. The improved performance due to used mineral oil, as boundary 
lubricant, is probably duo to some oxidation of tho oil and the formation 
of long chain acids in the oil. The advantages derived from adding small 
amounts of long chain acids to a mineral oil were recognized as early as 
1920 by Wells and Southcom.be [ 24 ] . A little later in 1921, Hardy 
and Doubleday [ 25,26] have reported tho effect of cha in length of 
additives in the lubricant. They observed a definite relationship between 
the coefficient of friction and the length of molecule ; longer the 
molecule lower the coefficient of friction. Beeok et.al [27 J , who 
used a very slow running four ball machine, showed that a coherent 
oil film, separating the surfaces, was promoted by the addition of a 
small proportion of oleic acid to white oil. 

Perhaps tho most telling support for Heeds observation is that 
provided by Fuks [28-29 ] . He performed a number of experiments with 

apparatus Very similar to that used by Heeds. The resemblances to Heeds* 
data are striking. It is of interest that an addition of o.l$ of stearic 
acid to a chromatographically separated naphthene-paraffin fraction oil 
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(MS-20) increased the residual film thickness. It is equally interesting 
to note that lower molecular weight hydrocarbons such as benzene, hexane, 
cyclohexane, isooctane and decane, exhibited no residual film. Fuks 
investigated many properties of the film and it was observed that film 
thickness increased with fatty acid concentration, with surface energy 
of the solid substrate and with increasing hydrocarbon chain length 
of both solvent and fatty acid. 

The works of Needs [ 22 ] and Fuks [ 28-29] have been a matter 
of great controversy and criticism for a long time by many workers, most 
recently by Hayward and Isdale [ 30 ] • It has been suggested that 
asperities on the solid surface might account for the residual film. Also 
dirt and undissolved material in the fluid have been invoked as 
explanations. According to Drauglis et.al [_3l] (who conducted a series 
of experiments with utmost care and observed residual films similar 
to Needs and Fuks), if the data of Fuks are all to be explained on 
the basis of dirt, then one must postulate dirt which changes its 
properties with temperature, length of additive molecule, concentration j 
of additive, but such a dirt would be a most unusual material. | 

Der jaguin L 32-33 J , who is regarded as the leading modern ; 

exponent of the view that molecular orientation in depth occurs near | 

f 

i 

solid boundaries and gives rise to significant rheological effects, 1 

j 

devised an experimental method in order to avoid any possible interference 
by dirt particles. Increase in viscosity was observed where traces of 1 

I 

1 

polar or surface active substances are present (including ordinary lubrical 

j 

oils) . ! 
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llore recently Cameron and his c (workers [ 34 - 37 ] have also 
studied the influence of surface active compounds in various lubrication 
systems. The lubricants were pure paraffins, mainly hexadecane, and the 
additives consisted of various long chain polar compounds. It was found 
that the oil film was markedly influenced by the additive. The scuffing 
load was found to increase with chain length of additive and solvent* 

Squeeze film studies similar to Fuks [ 28-291 were also made. The 

-4 

plates came to rest at a separation of about 2 x 10 mm, with all 
fluids, even with pure cetane. In addition, it was found that the 
increase in surface viscosity was greatest when additives and carriers were 
matched in terms of chain length and shape. The results presented 
in these studies are most easily explicable by postulating the presence 
of a surface film some 10' ” mm to 10 mm thick which has a much larger 
viscosity than the bulk of the fluid. 

Thus, although there may bo some controversy, as to the existence of 
rheological anomalies (like viscosity enhancement) in the neighbourhood 
of solid boundaries, the evidence on whole _seemfl positive. It is -evident 
that such zones of enhanced viscosity, giving rise to a residual film, 
would have important effects in the lubrication of practical machinery. 

Ordinary mineral lubricating oils contain various polar compounds, 
many of them surface active, in a non-polar hydrocarbon medium and 
thus belong to the class of liquids for which the evidence of the 
effect, viscosity enhancement , is strongest. The rheological anomalies 
could therefore have a considerable effect on the minimum film 
thickness, which in turn influences wear and fatigue failure of machine 


components 
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1.6 DEVELOPMENTS IF THEORETICAL PJIEOLO&Y 

A proper theoretical explanation of the residual film, reported 
by various ?rorkers, is clearly needed for a better understanding of 
thin film lubrication technology. An attempt in this direction is the 
work of Allen and Irauglis [ 9 ] 5 who developed a theory of boundary 
films based on the ordered liquid model, in which the wall forces, 
though not strong enough to bind more than a few molecular layer, induce 
long range order in the fluid. In this model, films on solid substrates 
are assumed to behave as liquids of special type, that is, they can flow 
as an ordinary liquid but have an internal order muoh like that in a 
crystal. This model fits qualitatively for all of the data of Fuks [28-29 j . 
A beginning of a quantitative theory based on molecular consideration 
is made in [ 9 ] a &d pursued farther in [38] • 

These detailed studies, although necessary for a fundamental 
understanding of physical and chemical phenomena involved in the formation 
of these films, are hardly justified for an engineering purpose. In 
such situations it may be sufficient to take into consideration the looal 
degree of freedom in some reasonable average sense. For example, the 
local rotation and spin of clusters and their stretch, in addition to 
classical degree of freedom, associated with a oontinuum may provide an 
adequate next step. 

Motivated by these considerations, attempts have been made by 
many researchersto extend the range of applicability of continuum mechanics 
(particularly in the study of fluid behaviour of rheo logically complex fluids 
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in situation where micromotion plays an important role) using the field 
theory approach to describe the macroscopic manifestation of microscopic 
events, i.e. microelement motions and deformations. 

The inadequacy of classical continuum approach to describe the 
mechanics of complex fluids has led to the development of theories of 
microcontinua in which continuous media are regarded as sets of structured 
particles possessing not only mass and velocity but also a substructure 
with which is associated a moment of inertia density and micro-deformation 
tensor. This extension of fluid mechanics required a reappraisal of 
classical concepts (e.g. the symmetry of stress tensor, the absence 
of couple stresses, etc.) in order to account for local structural 
aspects and micromotions. In fact, while many of the principles of 
classical continuum mechanics remain valid for this new class of fluids 
they had to be augmented with additional balance laws and constitutive 
relations. The presence of microscopic elements in a fluid gives rise 
not only to classical Cauchy stresses but also to couple stresses due to 
the microelement interactions. Further the orientation of the elements 
comprising the microcontinuum had to be considered as well as the 
production of internal angular momentum due to their rotational motions. 

No attempt is made in the present thesis to give an account of 
the historical development of the theory of microcontinua, as an excellant 
review has been given by Ariman, Turk and Sylvester [ 39 ] » and a reference 

may be made to that . 


The earliest formulation of a consistent general theory of fluid 
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microcontinua is attributed to Eringen [ 40-41 ] , in which the mechanics 
of fluids with deformable microelements are considered. The basis of 
Eringen’ s theory of "Simple Microfluids" was his work with Suhubi [42-45] * 

Although these papers were concerned with a nonlinear theory of microelastic 
solids, their treatment of motion, balance of moments, conservation of 
energy and entropy production is applicable to all continuous media 
where microelements are considered. 

In these papers [ 42-45] the authors develop a physical 
model in which each continuum particle is assigned a substructure, 
i.e. each material macrovolume element contains microvolume elements 

which can translate, rotate and deform independently of the motion of 

i 

the macrovolume ; however each deformation of the macrovolume element 
can be expected to produce a subsequent deformation of the macrovolume 
elements. Thus a mechanism is provided in the theory to treat materials, 
which are capable of supporting local stress moments and body moments, 
and in addition are influenced by the microelement spin inertia. 

The field equations (i.e. balance of momenta, conservation 
of mass and energy etc.) of Eringen* s microcontinuum mechanics have 
been developed by applying Cauchy* s law of motion and the classical 
conservation statements of mass and energy in the forms of volume 
and surface integral operators to each microelement contained in a 
material microvolume element. Simple statistical average of these 
global balance and conservation laws have been taken over the entire 
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macro volume element, resulting in the fundamental laws of microcontinua. 

The laws of classical continuum mechanics have "been augmented with 
additional equations, which account for conservation of microinertia 
moments and balance of first stress moment which arise due to consideration 
of microstructure. The key point to note in the Eringen’s theory of 
microcontinuum mechanics is the introduction of new kinematic variables, 
e.g. the gyration tensor and microinertia moment tensor, and the addition 
of the concepts of body moments, stress moments and microstress averages, 
to classical continuum mechanics. 

The mi cro fluid theory [ 40 ] is however too complicated and 
the underlying mathematical problem is not easily amenable to the solution 
of non-trivial problems in this field. This led Eringen [44] to 
postulate a subclass of these microfluids, called micropolar fluids, which 
would still exhibit the effects arising from particle micromotion but 
would produce considerable mathematical ease. Under these assumptions 
deformation of the fluid microelement is ignomed. Nevertheless, microrotatic 
effects are still present and surface and body oouples are permitted. In 
this theory two independent kinematic vector fields are introduced, the 
vector field representing the translation velocities of the fluid 
particles and the vector field representing the angular (or spin) velocities 
of the particle (i.e. microrotation vector) . 

Many theories similar to Eringen’s micropolar fluid theory have 
been advanced by various research workers in recent years from different 
view point. However the resulting field equations oan be reduced to 
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those of Eringen [44] (for a detailed comparison see rof. 39) • 

Physically, the theory of micropolar fluids may serve as a 
satisfactory model for description of the flow "behaviour of polymer 
fluids and fluid suspensions, where micromotion may play an important 
role* Another class of problems where micropolar fluid theory may 
provide an answer is' when the wave length of an acoustical disturbance 
becomes comparable to the mean grain-siae or distanoe in a body. In 
such situations diffraction patterns appear which can not be reconciled 
on the classical mechanical grounds « Under these circumstances the 
molecular or granular constituents of the medium are excited individually 
and the intrinsic motion of the material constituents must be taken 
into account. This situation could be thought of as one, representing 
the phenomena associated with thin films in proximity of a solid surface. 

1.7 MCROCONTINIJUM ANALYSIS OF FLUID PLOWS 

The theory of micropolar fluids has been applied extensively 
to describe the non-Newtonian behaviour of real fluids such as colloidal 
solutions, fluid suspensions and sometimes blood. A review of all such 
published work till 1974 was made by Ariman et.al, [45] » who listed 

over a hundred references. Since then soores of other papers have been 
published and no attempt is made to list them in a chronologioal manner. 

This section is restricted, basically to the discussion of those 
flow problems, which give evidence of rheological effects in the flow 
of micropolar fluids through narrow passages . 
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Increased viscosity due to the presence of microstructure in 
the fluid have been reported by Kline and Allen [ 46-48] ■ Cowin [ 49 ] 
and coworkers [ 50-51 1 » while studying the polar fluids [ 49 ] ( a fluid 

identical to Eringe^s micropolar fluid), have analyzed a special class 
of flows, called creeping flows. The most significant result of their 
work is that the polar fluid exhibits a resistance to motion, greater 
than or eq.ua 1 to that of the viscous fluid in the traditional flow- 
situations. Willson [ 52 ] and Peddieson and MciTitt [ 53 ] have also 
demonstrated increased effective viscosity in boundary layers in micropolar 
fluids . 

The origin of rheological abnormalities is then to be investigated 
on the basis of the theory of micro-mechanical motions, as there is a 
distinct possibility that micropolar fluid theory, which takes into 
account the intrinsic motion of the material constituents, may provide 
answer to some of the questions concerning these thin films, at least 
in a macroscopic sense. 

Motivated by theso considerations lubrication theory for 
micropolar fluids has been developed by many authors. The earliest 
work in this direction is by Allen and Kline [ 54 1 » who, using their 
own model, gave an approximate solution of the slider bearing problem . 
Similar theories were developed by Agrawal et.al. [ 55 ] » Balaram and 
Shastri [ 56 ] , Shukla and Isa [ 57 ] and Prakash and Sinha [ 58 ] 
using Eringen* s [ 44 ] model for micropolar fluids. 

Application of these theories has been made to various traditional 
lubrication problems. The slider bearing problems have been solved by 
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Datta [ 59 i » Maiti [ 60 3 , Shukla and Isa [ 57 J » Isa and 
Zaheeruddin L 61 ] , Agrawal [ 62 ] and Siriha and Singh l 65 J « 

Prakash and Sinha [ 58 , 64,65 3 , Mahanti [ 66 ] and Isa and 
Zaheeruddin [67] have solved the problem of one-dimensional 
journal hearings. The squeezing flow in parallel plates and curved 
plates have been solved by Agrawal et,al. [55] » Balaram [68] , 

Prakash and Sinha [ 69~7C>3 , Ramanaiah and Dubey [ 71 3 , Zaheeruddin 

and Isa [72 j and Tfyidon and Jaggi [ 75-74 j • Prakash and Sinha [75 3 
and Sinha [ 76 J have studied the dynamically loaded infinitely long 
journal bearings. Hydrostatic and externally pressurized bearings 
have been analyzed by Khader and Vachon [ 77 3 » Shukla and Isa [ 78 3 
and Isa and Zaheeruddin [ 79 3 « Prakash and Christensen [ 80 3 and 
Sinha [ 0l3 have analyzed the problem of rolling contact bearings. 

The elastohydrodynamic zone of roller bearing has been analyzed by 
Prakash and Christensen [823 , The inertia effects in squeeze 

bearings and thrust bearings have been studied by Mahanti and 
Ramanaiah [ 85 3 • 

The earlier theories were based upon the assumption that the 
Newtonian viscosity and micropolar viscosity coefficients are 
constant across the lubricant film. Recently Tipei [ 84 ] has 
given a more generalized form of Reynolds equation which assumes 
the viscosity coefficients to be function of all co-ordinates. 

However, Tipei [ 84 3 "too, has solved problem of short bearings 
under the assumption of constancy of viscosity coefficients. 
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The results obtained by using micropolar fluid as lubricant 
are encouraging and illustrate the capabilities of this theory in 
explaining the phenomenon of the enhancement of effective viscosity 
in thin films. However, the results in general are grossly 
exaggerated due to the inadmissible values of the governing parameters. 

Even after a decade this theory seems to be in a preliminary 
state. Scores of lubrication problems have yet to be solved from 
this view point. For example, little work has appeared on dynamically 
loaded bearings and porous bearings from the micro continuum view point • 
Moreover, the problems of finite dimensional bearings and bearings 
with rough surfaces from this view point are yet to be analyzed. 

Hence, in this thesis, an attempt is made to study these 
lubrication problems from the micro continuum point of view. 
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CHAPTER II 

BASIC EQUATIONS AND GENERALIZED REYNOLDS EQUATION 
FOR MICROPOLAR FLUIDS 

2 .1 INTR ODUCTION 

The theory of microfluids introduced, hy Eringen [ 40 ] deals 
with the class of fluids which respond to certain microscopic effects 
arising from the presence of microstructure and are influenced by spin 
inertia. This theory was later simplified to the theory of micropolar 
fluids [ 44 j in which the deformation of microelement is ignored. These 
fluids contain a suspension of rigid particles with individual motion 
which support stresses and body moments [ 40,44 J * This subclass is 
obtained from the general microfluid theory by imposing skew symmetric 
property of the gyration tensor, in addition to a condition of microisotropy. 

In view of the inadequacies of the classical Newtonian theory, 
lubrication theory for micropolar fluids has been developed by many authors 
[ 54*55*56,58 ] for two dimensional bearings, and by Shukla and Isa [ 57 J 
and Tipei [ 84 ] for three dimensional bearings , Apart from the analysis 
of Prakash and Sinha [ 58] , for the two dimensional case, no paper gives 

a detailed order of magnitude study for the three dimensional problem. 

Hence in this Chapter, a detailed order of magnitude study is made and a 
generalized type of Reynolds equation is obtained for the three dimensional 


case , 
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2.2 BALAHCE EQUATIONS FOR MICROPOLAR FLUIJDS 

The "balance equations governing the "behaviour of micropolar 
fluids have "been derived "by Eringen [ 44] and the results are simply- 
presented here. 

(i) Field Equations 

The field equations for micropolar fluids in vectorial form 
are [ 44 ] 

|f + V .(pv) = o (2.1) 

(X+ 2V) VV.V _ ± (2p+x) V * (Vx V) + XV*v~viI + pF = p — ■ (2.2) 

y . y . , y y . **>*->- _y. ->.)(. TY\) . 

(a+e+r) vv . v (vxv) +X 7 x V - 2x v + pL =p j — (2 . 3 ) 


The first equation is the principle of conservation of mass 
and the others are the conservation of linear and angular momentum, P 
is the mass density, V is the velocity vector, v is the microrotation 

y/ 

velocity vector, n is the thermodynamic pressure, F is the "body force 
per unit mass, L is the "body couple, 3 is the microinertia const ant, X,u 
are the viscosity coefficients of the classical fluid mechanics and X, a ,8,Y 
are the new viscosity coefficients for micropolar fluids . 


denotes material differentiation. 

**■ * 

For incompressible fluids V.V = 0, and n is replaced hy the 
hydrodynamio pressure p, to "be determined from the "boundary conditions. 
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The quantities X,u,x »«}&}! and j are constants and characterize 

the particular fluid under consideration. It may he emphasized that 

the constant X links together, through these field equations, the 

velocity vector Y, and the microrotation vector v. If mass, time 

and length are denoted hy M,L and T respectively, then X,y, and X are 

all of dimension I£L~ J 'T“\ while cx,g andy are all of dimension MIT" 1 

2 

and j of dimension L . 

(ii) Constitutive Equations 


Constitutive equations for the stress tensor t^ and the 
couple stress tensor ,are given in cartesian co-ordinates as 

- <- "* + + 2 X > + h,k> 


* ,k "" e k£r V 3p 


\ £ = av ^ + BV V 0 +Yv 0 


r,r kA 


k,£ £,k 


(2.4) 

(2.5) 


where E kjlr is an alternating tensor. An index followed hy a comma 
represents partial differentiation with respect to the space variable x^_< 

(iii) Boundary Conditions 


The boundary conditions are 

vOL^t) = Y b (2.6) 

v(Xg,t) = v B (2.7) 


where ^ is a point on the solid boundary having prescribed velocity Yg 
and prescribed microrotation velocity Vg .The zero spin or hyperstick 
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condition, Cowin [49] and Allen and Kline [ 54 ]» re pres en bo the 
extreme case in which the wall acts to completely inhibit microrotation. 
The zero spin boundary condition permits recovery of classical results 
by passing to the limit as a suitable measure of substructure size 
tends to zero. 

(iv) Three Dimensional Flow 


For subsequent study, attention is restricted to three 
dimensional steady motion of an incompressible micropolar fluid under 
the following assumptions 

(1) Body forces and body couples are negligible 

(2) all the characteristic coefficients are y-independent . 


Thus the velocity, the microrotation velocity and the pressure 
have the forms 


T = (u(x,y,z), v(x,y,z), w(x,y,z)) 
v = (v 1 (x,y,z), v 2 (x,y,z), v ( x ,y,z)) 


P = p(x,y,a) 


The balance equations (2.l) - (2.3) reduce to 


3u + 9 v + 3w _ 

3x 3y 3 Z 


|(2u + X )( ^ + 

3x 


„2 .2 
iJi + + x ( 

2 + 2 
3y 3z 


3v. 


37 
3 u 


3v 


3 P 


az ' a x 


f o u , „3u 3 un 

^ 3x 3y 3 z J 


( 2 . 8 ) 
(2.9) 
(2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 



*2 s 2 * * * * * a 2 3v l 

(s* 4X) ( M> + *4 + +*(rrk—£ --S 


dx By bz 


dz btJ By 


/ 3 V . 3 V , 8 V\ 

= o (U ■“ — + V + w — ) 

3 x B7 3 z 


1 (2U +x) £h + if! + 8 ±) + x (Jli . .! 1) . -ai 

2 3X 2 V V 8X 82 


/ 3w 3 w B¥ n 

- 0 (u — + ^ t — ) 

2 2 ^ 2. j 

, 3 V]_ 3 v i i- ,3% 3 V v 

Y ( — - + — - + — g) + x (gj- - *7£)~ 2X V 

3 x 3 y 3 z J 

3 V -, 3v 3v , 

3x 3y 3Z ; 

2 2 2 

3 v _ 3 v 0 3 v 0 

Y ( - 2 + 2 . + —) + Y ( SLH Qmm,) 2 X V 

T ' 2 + ~2 + 2' X Sz 3x ; “^ X 2 

3 x 3 y 3 z 

3v 3v . 3v . 

• ( 2 2 2 ' V 
= P 3 (u ~ — * + v — — + w ~~ — J 
3 x 3 y 3 z' 

2 2 2 
3 v, 3 v, 3 v, 

Y(_5 + _^ + _J) +X (-S._ffl) _2 X v ; 


3x 


3y 


3Z 


Pj( u 




3X By 


3v ■ 


DV -y 

4* y — 2 -f. -yyr — — 2) 

3x 9y 3Z y 


(2.13) 


(2 .14) 


(2.15) 


( 2 . 16 ) 


(2.17) 


2 *3 LUBRICATION PROBLEM-ASSUMPTIONS Ml BASIC EQUATIONS 

To state the lubrication assumptions a three dimensional slider 

type of hearing is shown in Pig. 2*1 as it is the simplest load carrying 

mechanism and an easier geometry to visualize the properties of more 

complicated films. U is the velocity of upper plate along 30-direction 

carrying a load W, h is the film thickness, agis the x-dimension and hgis 




Slider 

Oil film 

— » x 

Bearing 


Fig. 2.1 Slider bearing configuration 
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the z -dimens ion of the plate. The fluid is an incompressible micropolar 
fluid . 

The lubrication films to be analyzed are assumed to comply with 
usual assumptions |_ 85»86 j made in deriving the Reynolds equation 

1. Plow is laminar ; no vortex and no turbulenoe occur anywhere 
in the film. 

2. Body forces are neglected, e.g. the gravitational or magnetic 
field is assumed to be negligible. 

3. The film is thin enough as compared to the length and span 
of the bearing. This permits us to ignore the curvature of 
the fluid film and the rotational velocities may be replaced 
by translational velocities . 

4 • ho slip on the bearing surfaces • 

Apart from the above mentioned conventional lubrication 
assumptions, three more assumptions must be emphasized. Later on, in 
Chapters VI, VII and VIII, the first two will be violated. 

1. Bearing surfaces are smooth, i.e. there is no surface 
roughness * 

2 . Bearing surfaces are non-porous . 

3. The bearing surfaces are rigid, i.e. they are not deformed 
under pressure caused by fluid film. 

Equations (2,12) - (2.17) can be simplified with the help of 
above mentioned assumptions and following non-dimensional scheme. 
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3x 3y 3z^ 3y 3 a 


IP -n r - 5u , - 3ft , t JU 
— ■ f = Re {g,u + v — — + 6'n^“IT 

3x 3x 3y 3Z 


„ _2- a 2- p « 2- 0 3v n 9v, 

{ 2 Lz + Li t -« 2 - s -: }+ 2 h 2 {«, — i_ 6 , — 2 

1 35= 3y2 Z 352 2 35 1 32 


1 3p _ ri — 3v — ay ,f. — ay i 

£ = R e {6 u + v — — +S 2 w “ j 

1 3y 3x 3y ' 3z 


,r2 3 2 w 2 r c- *, 

1 a* 2 3 j 2 2 s5 2 1 -r. .=. 


2 3 v g 


ax 3y 


. i_ ni . Ee {5 a JS. + + 6 « -ii } 

5 3 3a 3x 3y 3z 


2 _ 2 - 2 - 

0 3 v, 5v, , 3 v, 

{6 2 i + t + 5, } + 


r 2 ! 2 


1 -=2 ay 2 2 2(l-N?) ‘ ” 2 


3x^ 3y 

A 2 


{-as - s 2 njL> 

ay 3? 


3 Vq 


^ — v -1 - Re {6,u 

(1 Jt) 1 1 


3X 


9v l 

+ t — + 5 2 w 

ay 


3v-, 


(2.18) 


(2.19) 


(2 . 20 ) 


(2 . 21 ) 


} ( 2 . 22 ) 
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2 2 - 2 

0 a v ' 3 vo 3 V 2_2 

J + + «2 ~ X + - iLi r«2 — - Si" 22 -} 

1 5,-2 2 ai2 2(lJS 2 ) 2 1 -- 


zy c 

' i^l 2 


3z 


3X 


- V 2 


( 1-N 2 ) 


v 2 = Re* {fi^u — — + 


3x 


{6 


^ t ix, 2 A 


A 2 

1 .* ' 3? 2 • “2 --- )+ - {S 


dVg, 

sy 

3Y 


3v, 

£ 

3Z 


+ « 2 w -} (2.25) 


3tt 


3x^ 


} 


3Z 2 2 ( 1-N 2 ) 1 3x 3 y 

3v -7 Bv-r 3 v . 


v - „ Re' {6 r -XI + T 

(l-J 2 ) 5 1 9 x 


— + (S 2 w 

3y 3 z 


> 

(2.24) 


where 


Re 


2 P hU 

2y+x 


Re' 


P>1Uh 

2 

4v* 


(2.25) 


The parameters N and k characterize the micropolarity of the 
fluid and distinguish a micropolar fluid from a Newtonian fluid. The 
parameter N, called coupling number, couples equations (2,2) and (2.5). 
For and these equations become uncoupled and (2.2) reduces 

to the usual Needier - Stokes equation. Thermo dynamio restrictions 

2 1 i 

require that o <_ N < ^ [ 84*87 J • 


The parameter k has the dimensions of length and can be thought 
of as a fluid property depending upon the size of the fluid molecule. 

In this sense the parameter L( = h/k) characterizes interaction between 
the geometry and the fluid. The case L -*■ « (in which the substructure 
size is negligible as compared to the bearing clearance) gives Newtonian 


results 



30 


It may be noted that Eringen [ 44 ] , Shukla and Isa [ 57 ] 
and Tipei [ 84 ] and other authors have used the parameter NL (where 
L = h/ £,) as one parameter, in discussing the micropolar fluids, hut 
this product characterizes a combination of effects, which are considered 
separately in this thesis. These parameters are discussed in more 
details in next Chapter. 

Re is the modified Reynolds number where u has been replaced 
by V + ^ X and this number is always smaller than the classical Reynolds 

x 

number. For lubrication problem Re <<c l, generally of the order of 10 ' . 

1 

Also it is reasonable to assume that Re << 1, since 3, is the square 

2 

of the length typical of the microstructure and comparable to l . If 
a requirement that microrotation velocity identically equal to the 
local fluid vorticity should be a possible solution of the balance 
equation, is imposed, it can be shown that 3 = % [ 47 ] • This requirement 
seems to be a natural one and has the desirable effect of reducing the 


number of constitutive coefficients. 

-3 

Further, 6^ and are of order 10 , 6 ^ is of order unity. 

Using these order of magnitude study and the above requirements, eqns . 
(2.19) to (2.24) reduce to 


— + 2N 2 
3y 2 




( 2 . 26 ) 


Pi = 
sy 


O 


(2.27) 
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^2 — 
o w 


1 

-aS- o 

C\J 

co 

ay 

' r 5 

?z 

2- 

+ jf/ 

SC . 

3w _ I^L 2 

s y 2 

2 (l-N 2 ) 


ay (i-k 2 ) 

. 2- 

iV 






ay 2 

1-ST 2 

* V 2 

- 0 

2 _ 

B 

P 2 ^ 

au , 

^L 2 „ 

ay 2 

2(l-U 2 ) 

ay 

(UN 2 ) 


(2 . 28 ) 


(2.29) 


(2.30) 


(2 .31) 


From equation (2.27), P = p(x,z), and when (2 ,30) is solved with 
zero boundary conditions [ 49* 54-] » v 2 == 0 * 

(i) Velocity Distribution 


To solve for velocity components and microrotation velocity 
components, we revert equations (2.26) - (2 .31) to corresponding 
dimensional form . 


The equations are 

. r.2 9v, 9p 

( p+ l x ) + x __1 ^ > = 0 

" ay ay a* 


, -2 av. 

(„+ Iv ') LZ _y 1 _ _ 0 

ay2 x ay az 


2 ' 
a 2 vi 


J 1 a w 

Y ?r - 2Xv. +x — — - o 

2 1 3 y 

ay 


(2.32) 

(2.33) 

(2.34) 
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9 2 V3 

Y 1 - 2 xv 

s/ 5 



o 


( 2 . 35 ) 


The four differential equations (2,52 to 2.35) form two 
systems of simultaneous equations. One consists of (2.32) and (2.35) 
and the other of (2.33) and (2.34). The solution of first system will 
yield expressions for u and v . The second system will yield expressions 
for w and v^ , it is seen that the second system "becomes identical to 
the first system if we replace "by -v^. Hence it is necessary to 
solve only the first system for u and and the expressions for w 
and v x can "be written in a similar way. 

The boundary conditions for velocity and microrotation 
velocity are 

n = U-j_i> w - ^12* v i ~ °» v 3 ~ ° ,at y = 0 (2.36) 

U = U 21’ W = U 22’ V 1 “ °* v 3 “ °» at y ~ h (2.37) 

The solutions are 

2 -^2 

u = ~ (“3^ 2 + A lx y) - — tA 21 sinh my + A^ cosh my} + A (2.38 ) 

Vj = - 2y — 7 + ^ll) + A 21 cosh + A 31 siah “y ( 2 *39) 

T -£$ 2 + h2 y) - f- { A 22 sinl1 + A J2 oosh ”y }+ A 42 

(2.40) 

Vj - - -f - ££•( -jf y + A 12 ) + A 22 cosh ay + A^ 2 sinh my} (2.41) 
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where A^, i = 1,2, 3,4 i j a 1,2 | are constants and given by 


A, . » [ (lT 0 . - U .) sinh mh th sinh mh + ~(l - cosh mh) } 3 /A 

Xj 13 ^ 3X. ™ 5 

J 


A . 

a = -la 

2j 2y 

1,. - i[ (0, .JJ .) b - - °°MH . #i) + £ £B_ { £ ( oosh mh -X) 
3j y 2o „ 2n 3x. 2 

3 


(2 .42) 
(2 43) 


2 

+ (h - ~ sinh mh) } 3 /A,_ 


A 4j 


2N 2 

TJ_ . + ~ A, . 
10 m 33 


2H 2 


5 "y 


{ sinh mh - — £ (cosh mh - l)> 


(2.44) 

(2 .45) 

( 2 . 46 ) 


in which 


x. 


N 


^ = x, %2 - z, bi = “ 


(2.47) 


(ii) Flow Flux 

The volume flow rates and q. z , along x and z axes, are 


given by 


Sc* l udy 


(2 48) 


and 


% “ / w dy 
2 -'o 


(2.49) 
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or 


a (V+V,) 

X 2 v ll 21' y 3X 


(2.50) 


and. 


\ ~ 2 " U 12 +TI 22 


) _ h^f(hjl,h) jP 

y 3 2 


(2.51) 


where 


2 

r* /rr n , \ 1 , jHjl ,, ITh 

f ( N » £,h ^ “ 12 2 “ 2h coth 2?. 

h 


(2.52) 


It should be noted that for convenience in Chapters VI and VII the flow 
fluxes are defined as follows 


„ £(„ + o ) _ £PM.h) 21 

x 2^ U 11 U 2l' y 3x 


« ktjj +u ) _ Klal A ), -11 

OV U OT U 00' f 11 -Mr 


bs 2 V 21 22 


3Z 


(2.53) 

(2.54) 


where 


„/, T „ , \ h^, , .2 N £ h 2 ,, Kh 

f(h,£,h) = ^- 2 + h A - -g — ooth 


(2.55) 


(iii) Generalized Reynolds Equation 

Integrating the equation of continuity (2 .ll) across 
the film and using (2.50) and (2,51) (or eqn s«2»53 and 2.54) 
together with (2.36) and (2.37) the generalized Reynolds equation 
for micropolar fluid in three dimensions is obtained as 


3 <«*> + lbr K 


3x 


-12l + p -H + v -V 

21 9x 22 32 oh 


(2.56) 
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where and V are the normal velocity components of the surfaces 
at y = h and y = o, respectively. 

Similar equation was obtained by Shukla and Isa. [ 57 ] 

and Tipei [ 84 J under the assumption of constancy of characteristic 
coefficients across the film. However the parameters iT and ~ 
are different from the parameters used by [ 57»84 J • 

For infinitely long bearings, along z-axis, the two 
dimensional Reynolds equation is given by [ 58 J • 


(iv) Frictional Torques 


The shear stresses along x-direction and z-direction, 
from equation (2»4)» are 

- (w+ f X) ~g + XV, 


X 


(2.57) 


and 


% ■ (“ + 2 1 X) ^ ^ 


(2.58) 


and on the surfaces y = o s,nd y = h these assume the form 


( x ) , 

x'o,h 


- h n ^21^11^ 

+ 2 3X + hg(H,£ ,h) 


(2.59) 


(t ) 

' z'o f h 


- h 3P V ^ J 22 JJ 12^ 

+ 2 dz hg(N,£ ,h)" 


( 2 . 60 ) 


where 


g(h, Jt,h) 


i _ m tanh m 

1 h * m 2£ 


( 2 . 61 ) 
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The frictional drags are obtained by integrating the 
shearing stresses over the bearing surfaces. 

It should be noted that for convenience in Chapters VI 
and Til the shear stresses are defined as follows 


/ \ - h u ^21~ ,u n- > 

' x'o,h 2 3x g(N,£,h) 


, . - h 9 P y(U 22~V 

^Vo,h + 2 p + g'CMjh) 


( 2 . 62 ) 

(2 .63) 


g(N,£,h) - h - 2N£ tanh 


(2.64) 


where 
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CHAPTER III 

TWO 4FD THREE DILEEL7ST OI'TAL IW-.CTCT.IC S0T3EZE FILMS 
3.1 INTRODTJCTIOH 

In recent years many workers have analysed the problems of 
micropolar fluid lubricated non - cyclic squeeze films [ 55»68-71,83 J 
and cyclic squeeze films [ 75-76 ] . The problem of porous non-cyclic 
squeeze films have also been solved using micropolar fluid as lubricant 
[ 72-74 ] • The prominent feature of using micropolar fluid as 

lubricant was an increase in load carrying capacity and squeeze time. 

Though some of the geometries, such as circular plates, 
considered were finite dimensional, but due to non-variation of 
pressure along circumferential direction, the problems solved v,rere 
essentially two dimensional. It is thus imperative to study the 
effects of micropolarity on three dimensional squeeze films problans. 

In this Chapter, an application of micropolar fluid theory 
is made to study some hitherto untouched two dimensional important 
problems and to some three dimensional configurations. 

The problem considered is that of the so called T non-cyclic 
squeeze films’ . The squeeze velocity V is imparted by a load W, which 
will be considered constant in magnitude and direction and symmetrically 
placed, with respect to the boundaries of the system. Various geometries 
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are considered , namely , spherical hearings, slider and rectangular 
plates, conical hearings and some miscellaneous configurations. The 
approach is similar to that of Pinkus and Sternlicht [85] * Detailed 
numerical results are obtained for the spherical hearings, whereas 
for the three dimensional case only the dimensionless response time 
is obtained. 


3.2 SPHERICAL BEARINGS 


The film geometry for squeeze film in spherical hearing is 
shown in Fig. 3.1(a) • Here we shall consider the case of hemispherical 
seat . 

From the geometry, given by Fig .5.1(a), the film thickness 
is given by 


h <=* r - R - a^cos e (3*l) 

which can he rewritten as 

h *= c(l -e cose) (3.2) 


where 


o = r -K, s - 

For a sphere the amount of fluid passing through a conical element 
of surface is given by 


(3.3) 


h 

q « / (2 ttR sine)u dy 

r\ 


(3.4) 


where u is the velocity of the fluid. For zero tangential velocities 
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the expression for u, equation (2 .58), reduces to 


u 


„ _1_ z 2 A oosh my 1 \ h dp 
nE d9^ 2 ” m s inh mh “ 2 pR A0 


IT 

7 - ~ * 

m 


x { sirih my - ^ } 

sirih mh 


(3.5) 


in which x is replaced by Re . 

Substitution of u, equation (3.5) in equation (3*4-)* yields 


q 0 = 2 tt R sin 0 ~ f(l,*,h)} 


(3.6) 


where m,N and l are given hy eqn. (2.18) and f (lT,£ ,h) is defined 
by eqn. ( 2 . 52 ) . 

The displacement of the fluid at any point in the seat is [ 88 ] , 

q - *VR 2 sin 2 e (3.7) 

0 

where V is the relative normal velocity of the surfaces, responsible 
for squeeze film action and is given by 
de 


7 - c 


dt 


(3.8) 


Equating (3.6) and (5.7) 

2 

dp s uVR sin 8 
de 2h 5 f(R,S',h) 

The pressure distribution with the boundary conditions 


(3.9) 


p( +2) “ 0 


( 3 . 10 ) 


becomes 
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ir 

J2. 7 


p(s) - ^ , OO < 

9 h 3 f(lt,S ,h) 


IT 
1. 2 


(3.11) 


V 

since the pressure distribution is symmetric for the regions ~ 9 <_o 

and o ^_0 <_~2 • 


Non-dinencionalizing, using 
H = — = (l- e cos e) , i = - , P = 


VlR 


2o 


(3.12) 


the non-dimensional pressure is 

7T / 2 


p(e) = / 


sin© d0 


e h 5 p(n,l,h) 


(3.13) 


where 


F(N,L,H) - -■ 


1_ 1 N 

12 (LH) 2 “ 2HI 


coth 


HLH 


(3.14) 


The load capacity for a hemisphere is 


2 f 7 ’’/ 2 Q TryTR 4 / 1 / 2 sin 3 9 d! 

= 2v R ; p sm 0 cos 0 de - -~x — ■ J 

n 21 


0 h 3 f(u,&,h) 


(3.15) 


The non-dimensional load is 


W 


(=^) 

2c 9 


= j- 17 / 2 sin 5 6.d0 


0 H 5 P(H,L,H) 


(3.16) 


The time of approach, t, from an initial film thickness 
h x to a final film thickness h Q (h 1 > h Q ) in terms of corresponding 
values of eccentricity ratio can be evaluated using (3.15) and (3.8). 
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Thus 

e i tr/ p 3 
± j ' d sin 9 d8 de 

0 0 ' b 3 f(lT,Ji,h) 

The non-dimensional time of approach is 

ip _ t _ j sin 3 9 d9 de 

(JLESi.) e o 0 h 3 f(n,l,h) 

2o 2 W 

3.3. SLIDERS A HD RECTANGULAR PLATES 

The rectangular plate poses a three dimensional problem, 
and we must return to the basic differential equation of a viscous 
flow of micropolar fluid for a solution. 

The differential equation governing the pressure distribution 
in three-dimensions, equation (2.56), for h = constant, is 

iViit —bi ( 3.i9) 

3x 3Z h' ) f(N/,h) 

Equation (5*19) is Poisson's equation, the solution of 

which is assumed of the form 

00 

p(x,z) = l (A cos -2^-2) 2(a) (3.20) 

n=l,3» 5 n a 2 

where Z(z) is the function of z alone and A^ is a constant to be 
determined . is the x-dimension of plate. Substituting (3 .20) into 
(3.19) and using the relation 


t = 


jrUclT j 

e 


2W 


(3.17) 


( 3 . 18 ) 
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00 (n-l) 

\ r (-1) 2 

n=l,3>5 


1 nux 

— cos ~~~ 
n a 0 


we have 


l (— ) 2 A oos^iSzW + 
n=l,3 f 5 v a 2 ' n a 2 


I A uttx d Z (z} 

L JL cos 

n«l,3j5 11 


a 2 dz 2 


_ 1 


PV 


" ^5 


i . . (-D 


(n=ll 


h 5 f(N,%h) n=1 »5»5 


1 iittx 

— COS — 

v. a o 


which is satisfied if 


A. 


A£M-(^) 2 A n z W -^ 


(n-ll 
2 


n dz 2 


A particular solution is given by ^ 

mz _ , njrz ^ a 2 V — .Irll.. 

Z (”) = B m sinl1 AI n *, ' A 3 f (H,t ,h) n\ 

Putting (3.3$ into (3.20) and using the boundary conditions 

\ 

p = o at s *» + 2 

where b 2 is the z-dimemion of the plate and hy the requirement 

symmetry in z-direction, the pressure distribution is 

JScl 


nrrz 

a 


+4U a 2 V “ (-1) 2 [ C0Sl:1 

-1,3.5 u 3 


- 1 


cosh 


2a, 


(3.21) 

(3.22) 

(3.23) 

il 

(3.24) 

(3.25) 
of 

] (3.26) 
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The total load is given by 


¥ = + 


16 u a'Y 




* a nrrb ? b ? 

l \ -z ^ "T* tarih "~2aT “ rr 
n=l,5»5 nrir 2 2n 


) 


(3.27) 


dh 

Using Y = - ^ , the time of approach from an initial film thickness 


h^ to a final film thickness h Q (h^> h Q ) is given by 


h. 


dh 


16 v &L « b a 0 nv b xx i 

—7 I ( ~T “ T* tanh ’~~2sh / T") 

irTT n=l,5,5 2n 4 r?ir 2 h Q h 5 f(U,£,h) 


(3.28) 


Yon-dimensionlizing, using 

h 


h 

h. 


H - ¥ , H - 


V 1 " 


T 


( 


16 vi a' 
TT^Wh^ 


3 - I ; h_ h_ tanh 
11 - 1 - 1 ’ 3 * 5 ' 2n 4 ' A 


nvb , 
2a_ 


(3.29) 


the non-dimensional time of approach is 

T = J 1 (3.30) 

H q ErtKl,L,H) 

3.4 ELLIPTICAL PLATE5 

For an elliptic plate 



(3.31) 
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the pressure distribution with the help of governing differential 
equation (3.19) can he written as 


2 2 

p(x,z) = C n (—■ + —- l), C is a constant 
1 "1 


(3.32) 


where the pressure boundary condition is p = o on the boundary of 
the ellipse given by equation (3.3l). 

Substituting equation (3.22) into (3.19)»we get the 


pressure distribution 


p(x,z) 


2 2 
WA i B i 


2 2 

(Vi- 1 ) 


2 (A^+B?) h 5 f (N/ ,h) A 1 B 1 


(3.33) 


The load capacity and time of approach from an initial film 
thickness h^ to a final thickness h Q are given by 


W = + 


a.3^3 

A 1 B 1 


4h 5 f(H,^,H) ( A i +B i) 


(3 .54) 


* PA^bJ h 

1 11 f dh 

4W (A^+B^) h Q h 5 f(N,&,h) 


(3.35) 


Circular plates are special case of elliptioal plates and the bearing 
characteristics can be obtained directly from above expression by 
putting A 1 =B 1 =r, Prakash and S inha [69] • 


3.5 MISCELLANEOUS CONF I&UBATIQNS 


By the methods similar to those described above, results for 
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other shapes, of falling plates, can “be extracted. Some of these 
are given hy 


(i) Circular rings of inner radius r^ and outer radius r 2 

2 2 X 2 
r^ 2 ” r 

u 2 — i»L r 


F 


ir yy 


r a c\ a 

r „4 _4 r r 2 " r l^ -i 

3 , r 2 “ r l” r 

8h 3 f( K ,£,h) log fe 


(3.36) 


putting 


_ & 


V = - 


dt 


t = 


Try - /i /i ( r o~ r -t) 

q¥ , ^ 

log(^) 

1 


2 2x2 h 


* / 1 


dh 


h h 3 f (F,£ ,h) 


(3.37) 


(ii) Circular sector of angle g ^ and hounded hy radii r^ and r 2 


¥ 


8V6JT 


Tr 2 h 5 f(¥,Ji,h) 


4 4 

2 1 


l [ 

1 > 3 ’ 5 4 n 2 ( 2 -* i -) 2 


2~rr 


1 


U-rfM 2 ] 2 

1 


o+S3L 

2.nfr 2n f 

6 i 

r 2 1 - r x 1 


( 3 . 38 ) 


and 


8 3 x y 


t ~ 


I 

n=l,3»5 


4 4 

T—T 
2 1 


,nrr 


c- > 7> 

r„ e 


2 P 1 " r l P 1 


2 + 


tlTT 


. 2 / n^ \ 2 

4n ( 2 -g — ) 


2n-rr 

- 


2nir 

r P 1 
r l 


2tt 


n^U-tf 21 ) 2 ) 


dh 


I 

h h 3 f (H,£ ,h) 


X 


(3.39) 
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Letting r -»• o in equations (3.38) and (3*39) the characteristics for 
full sector can he obtained. 

(iii) A circle of radius r T , a -portion of 'which has been cut out 
by another full circle of radius r^ (fig. 3.1(b)) 


W 


MY 


r / 0 4 , 2 2 An -1 
L (2rIT-4r r -r ') cos 


8h 5 f(N,£,h) 


-l^ r2 — 2 


2r„ 


4 * 


+ r (r H * 2 ) (4*H> 1/2 ^ 


(5.40) 


and 


t - + 5T [ (2rJ-4r^r^) cos' 1 




h h 5 f(lT,^,h) 


Letting r 2 ->■ o , the case of full circle is obtained, Prakash 
and Sinha [ 69 ] * 

(t- v ) Full Conical Seats 

In case of complete cone (Fig. 3 .l(c)) » 


W 


tt\i R^V 

* Jt- 


8l4f(lT,i,h) sirrot^ 


(3.41) 


(3.42) 


where a-^ is the semi-vertical angle and is the radius of base 
and 


>4 h 


t = + 


I 1 


dh 


8Wsin 4 a 1 h Q h 5 f(N,£,h) 


(343) 
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Fig 


3 1(c) Full conical seat conf i g uration 
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(v) Truncated Conical Seats 


For truncated cone (Fig. 3»l(d)), whose smaller radius 
is Rg and larger radius is R^, 


W 


TTUY 


8h :J f(lT,ii,h) sin 4 ^ 


n 0 0 

f r 4 r 4 -i R r R 2 ) n 
L \~X-2~ R _ J 

!og(~) 


and 


t = 


r,4 „4 C 4^ )2 ■, rh ah 




ewsin^'” 1 " 2 ^ h 0 h 3 f(n,Jt,h) 

R 2 

The dimensionless response time in each of the above 
cases is governed by integral 


/ 


dS 


H - h 5 f(n,l,h) 


This integral is therefore integrated numer ioally for 
rectangular plates and represented graphically. This integral 
is similar to the one obtained by Prakash and Sinha [ 69 ] • 


3.6 RESULTS MR DISCUSSION 

(i) Dimensionless Parameters 

The micropolar fluid is characterized by three physical 
constants, V ,x and Y as compared to a Newtonian fluid which only 


(3.44) 


(345) 


(3.46) 



J I 



Fig 3 1(d) Truncated 


conical seat configuration 



52 


has one constant y, the Newtonian viscosity. This arises as a consequence of 

resistance to linear motion, v/hereas the parameter X , which has 

the dimensions of viscosity may he visualized, to be a consequence 

of resistance to the micromotion of the fluid molecule. The third 

2 • 

coefficient X, which has the dimension L L U J » gives rise to 

a very important parameter, when compared with the Newtonian viscosity y, 

y 1/2 

i.e. the parameter £ = (^j~) ' . This £ has the dimensions of length 

and may be thought of as a material parameter, called the material 
characteristic length, depending upon the molecular size of the fluid 
or the chain length of the rigid particle additive . 


Thus these three physical constants,y ,X , and Y give 
rise to two~dimensionless parameters defined as 


N - (- 


2h+x 


) 1//2 and 


L =f 


l_Nl/2 




( 5 . 47 ) 


called the coupling number and the length ratio. These parameters 
are not encountered in the usual Newtonian theory, and are likely 
to be of great value in discussing the effects arising due to the presence 
of rigid particle additives in the lubricant. 


The coupling number N, characterizes the ooupling of the 
linear and rotational motions arising due to the micro-motions of the 
fluid molecules or the additive molecules. Thus the number N signifies 
a coupling between Newtonian and rotational viscosity. Since the 
apparent or equivalent dynamic viscosity u is a derived magnitude 
including X and another independent constant y, v , oalled dynamic 
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viscosity, p = u v + — x [44] • Thus the expression for IT becomes 

IT 2 - —* 

2(W V +X) 

Thermodynamic restrictions require that , X 21° [ 40j44] , 

it is thus evident that o f_ IT 2 < ^ . 


(3.48) 


The second parameter L, defined as the ratio of the clearance 
width to the characteristic material length & , characterizes interaction 
between the geometry and the fluid and may be considered to be a property 
depending upon the size of lubricant molecule, say, the chain length 
of the molecule. 


It is expected that the effects of micropolarity will be of 
importance either when the characteristic material length is large, 
corresponding to a large size of lubricant molecule or when the 
clearance width is small. The second case is of importance in 
lubrication theory where c is usually small and hence comparable with & . 
Recalling that L = c/l and keeping c fixed an increase in a means at 
decrease in L, so lower the value of L, larger the chain length 
of molecule and higher should be the effects due to micropolarity. 

In fact it is seen that as L ■* o, the expressions for velocity and 
other flow characteristics reduce to their equivalents in the Newtonian 
theory, with U everywhere replaced by (u + ^x) , Thus indicating an 
increased effective viscosity. This observation supports the experimental 
evidence of increased effective viscosity t 8,22,28—36 ] in the presence 
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of additives having large chain length and as well as in the case of 
thin films. Since a small value of L corresponds also to a small 
clearance width provided £ is kept as constant . A large value of L, 
on the other hand, signifies a large clearance width, i.e. a case of 
hydrodynamic lubrication or a small material characteristic length 
and the effects of the micropolarity are not likely to he significant. 
It is seen that as L tends to ® , £ would tend to zero, and the 
individuality of the microelement is lost and once again the Newtonian 
results are obtained. 

(ii) Bearing Characterist ics 

Out of the various cases that have been treated in this 
Chapter, the numerical results are obtained in detail for the case 
of spherical bearings (Fig. 3.2 to 3«7) whereas for the three- 
dimensional rectangular plates only the response time is plotted 
(Figs . 3 *8 and 3 *9 ) • 

Non-dimensional pressure distribution curves are shown 
in Fig. 3.2 (N =» 0.7 and various values of L) and Fig.3«3(L = o.o 
and various values of N) • Since the pressure curve is symmetric 

r TT v 

about 9=o, only half of the pressure curve is shown (o <_ 6 <_-^j . 

L ■+ “(Fig. 3.2) and N o (Fig, 5.3) show the Newtonian pressure 
curves. As the value of L is decreased, which corresponds to a small 
clearance or a large size of additive molecule, the pressure begins 
to rise (Fig, 3.2) and for L = 0.0, the maximum pressure 

is obtained , As the coupling number increases 
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(Fig. 3.3) the additives play a significant role. The pressure in all 
cases is higher than the Newtonian pressure. 

Fig. 3.4 and Fig. 3*5 show the graphs for dimensionless load 
capacity as a function of micropolar parameters N and L. The load 
also is seen to increase as the value of IT increases and L decreases. 

The Newtonian curves are straight lines as the load remains unaffected 
by the variation in N and L in this case. 

Figs. 3.6 to 3*9 are the curves for response time for various 
values of N and L. Response time is actually the most important 
squeeze film characteristic. This gives us the time that will elapse 
for an oil film to be reduced to some minimum value, a value which 
may have been established as the thinnest permissible film for a 
particular bearing. The time required to reduce the film thickness 
can then be compared with the length of the time the load is actually 
applied. 

Fig. 3*6 shows the time-height relationship for a fixed 
value of N(N = 0 . 7 ) and various values of L. For large values of H, 
i.e* for large film thickness, the effect of additive on response 
time is not significant. But as the surfaces further approach each 
other, the value of H decreases and a marked discrepancy is observed 
between the Newtonian and micropolar curves. This discrepancy is conspicuous 
for H < 0.3 and for low values of L, A similar result was obtained by 
Needs £ 22 ] in his experiments on squeeze films. This discrepancy 
indicates that for low values of L and H, the time of approach is 
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increased, which, in turn indicates the enhancement of effective 
viscosity, which is a confirmed experimental fact './hen additives 
are used in lubricants. Similar interpretations can be made 
regarding fig. 5.7. 

Pigs. 3*8 and 3»9 depict the response time curves for 
various values of L and N for rectangular plate squeeze film 
(three-dimensional problem) . The qualitative effects of micropolarity 
on response time is same as that in two-dimensional problems which 
can be easily ascertained from these graphs. 

It can thus be concluded that the theoretical effects of 
additives on three-dimensional lubrication problem are almost 
identical to those in two dimensions, at least qualitatively. 
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Fjg. 3*2 Non dimensional pressure distribution for 
various values of L at N = 07 and G-0-7 
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|.3‘4 Non - di mcnsior-al load vs. L for various value 
at 6=0-7 


of N 
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Fig. 37 Non-dimensional height vs. non-dimensional 
time of approach for various values of N 

at L -2*0 


UJ ; 



time of approach for various values of L 
at N = 0 7 (Rectangular plate) 




Fig. 3-9 Non-dimensional height vs non-dimensional 
time of approach for various values of N 
at L = 2*0 (Rectangular plate) 
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CHAPTER IV 

DYNAMICALLY LOADED SHORT JOURNAL BEARINGS 

4.1 INTRODUCTION 

Dynamically loaded bearings in which the load alternates 
or rotates has been subject of much interest, conjecture and research. 
Dynamic situations where both wedge and squeeze films are active, arise 
in the bearings of almost all reciprocating machines, especially in 
steam engines, automobile and aircraft piston engines and various 
machine tools. Even the human joints can be approximated by 
dynamically loaded bearings* 

The load-carrying phenomenon in such bearing arises due to 
the fact that a viscous lubricant cannot be instantaneously squeezed 
from between the two approaching surfaces. Because of the resistance 
of the lubricant to extrusion, pressure is built up and the load 
is supported by the oil film. When the load is relieved or reversed, 
the oil film thickness can recover its thickness before the next 
cycle, if the bearing has been designed to permit this build up. 

This phenomenon is found in reciprocating machines where the 
bearings are subjected to fluctuating, impulsive or dynamic loads 
and therefore present problems not met in the normal operation of 
steady load bearings* In general, lubricating oils in the bearing 
subjected to reciprocating loads are highly vulnerable to contamination, 
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Toy both dirt and metal particles'. The lubricant may then behave as a 
fluid suspension. Moreover, the use of additives is a common practice 
not only for the enhancement of viscosity and a consequent increase in 
the load carrying capacity, hut also to increase the working range of 
temperature. These additives, through chemical reactions, with the 
metal surfaces form solid surface films and thus protect the hearing 
surfaces from damage and wear. 

Little work has been reported [75*76 ] on the lubrication 
of the dynamically loaded bearings from this view point. It is thus 
felt that the application of the concept of micropolar fluid lubrication 
to bearings operating under dynamic conditions would lead to a better 
understanding of bearing operation and could yield results in conformity 
more to practice than those obtained by the Newtonian theory. The 

i 

application of the present work may also be found in the lubrication 
of human joints, which can be approximated by dynamically loaded 
bearings. These joints are lubricated by synovial fluid, which has 
long ohain hyaluronic acid molecules and can thus be thought of as a 
micropolar fluid. 

In earlier works, by Prakash and Sinha [75] a nd Sinha [ 76 ] , 
the problems solved were restricted to infinitely long journal bearings. 
However, for practice, the short bearings are better approximation to 
finite journal bearings [ 89 ] ♦ Hence, in this Chapter, we solve the 
problem of dynamically loaded short journal bearings lubricated with 
micropolar fluid. 
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We start with the general case of dynamic loading in which the 
load and speed of rotation vary hoth in magnitude and in direction with 
time, Various hearing characteristics for an infinitely short journal 
hearings are obtained in the form of integrals. The equation for 
determining the journal locus is rather complicated and therefore only 
a particular case is considered in detail, name ly pure squeeze films 
with no journal rotation. The case of a fluctuating (sinusoidal) load 
without rotation is considered. Under this loading condition, the 
journal may he thought to act like a piston in hydraulic damping 
mechanism. Physically this implies that, as the load is applied 
in one direction, the fluid is squeezed out on one side and sucked in 
on the other. Before the piston comes in contact with the hearing, 
the load reverses its direction and the oil on the other side begins 
to take the load. Thus, assuming a continuous fluid film to exist 
from o to 2tr the hearing response to sinusoidal loading is studied 
and the micropolar effects are elaborated graphically. 

4*2 THE PROBLEM AND MAIYSIS 

The problem considered is that of dynamically loaded journal 

hearings (Fig.4»l) in which the load is a function of time. The journal 

is in combined tangential motion U and normal motion V , where U and 

V are hoth functions of time, 
o 

Equation (2 ,56) , governing the pressure distribution, for 
the type of flow mentioned above, assumes the form 
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3 


j r (jV, Af h)^> + 


{ ff(Hn,h) 

BZ y v ’ 1 J BZ 


3P 


+ h -H + v 

2 3x 2 3x o 


( 4 . 1 ) 


where N and £ are given by eqn. (2.18) and f(R, £,h) is defined by 
eqn. (2 .52) , 

The first term in the r.h.s. of eqn. ( 4 «l) represents the action of 
the journal rotating with a velocity U over a wedge shaped fluid film 
given "by h(x) , the second term implies a variation of tangential 
velocity along the hearing surface and the last term is the expression 
for velocity of the shaft centre and is responsible for squeeze film 
action. 


Assuming the bearing to be fixed, the shaft centre will 

have instantaneous radial and tangential velocities. Referring 

d e d (<f> +$) 

to the vertical line ip - o, those will be c and e u 


dt 


dt 


Thus U and V Q are given by 


U = Rrn + c ~~ sin 8 - c e 


cos 6 


( 4 . 2 ) 


and 

V = O'jjT" cos ® + 0 e s ^ n 6 ( 4 . 3 ) 

o dt dt 

Substituting for II and V q in equation (4.1) and considering 
that, L 85 i » x=R 9,|-«2 and ~ cos6« 2 and 
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2 re ili±il sin8 , || 

dt dt d6 


(4.4) 


one obtains, the Reynolds equation for the general case of dynamic 
loading for micropolar fluids, 

1_ 9 r h 3 f(ll,& t h) 3 P , JL_ { h"f(R,£,h) 3P , 

O rv /-\ <■ - , A A -> A r* .. A rj> -* 


R‘ 


,2 80 


80 


8 z 


3z 


/ 0 d4 \ dh, , n d £ Q 

i “ 2 dT> aT } + 20 dT 008 e - 


(4.5) 


For short journal hearings, eqn. (4*5) reduoes to 


_ _ — 
dz 2bf’f (lT,&,h) 


y { (w- 2 <*> t - 2 ~-) H + 2o,~- cos 9} (4.6) 


f dif; N 

where as is the angular velocity of the journal and = -J£“) is the 
frequency of rotation of applied load. 

Integrating eqn. (4*6) twice, using the boundary 
conditions 

L 


p( 1 i~) - 


(4.7) 


the pressure distribution is obtain ed as 


p(6,z) - 


4h 5 f(N,&,h) 


{ (u-2 w T - 2 ^-) •— + 2c cos 0 } (z 2 -7^) 


dt 


(4.8) 
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The load component, , normal to the line of centers is 


given "by 


^o 

" 2 " 2 11 

¥ = ¥ sin 4 = 2 / / p(0,z) sin0R d0 dz 

2 00 -' 


( 4 . 9 ) 


Substituting the value of p and integrating, we get 


W = ¥ sin 4^ 

7 T 


PeCl^R A A , 2 " - 2 fi _ fl 

0 / 0 0 d<£ Nf sxn 0. d9 

-(to - 2 03 T - 2 


24 


0 h 5 f(V,h) 


(4 . 10 ) 


The load component, W t along the line of centers, 


is given hy 

ib 

2 2TT 

¥ = ff cos tj> = -2 J f p(e,z) cose R d6 dz 
0 0 0 


(4.11) 


which after simplification "becomes 

oy A 2* 2 

rpr __ 0 d£ f GOS V CLP 

\'q no if J ^ 


12 dt 0 h 5 f(N,£,h) 


(4 .12) 


Equations ( 4 .10) and ( 4 .12) give a relation "between 
the applied load and the motion of shaft center. In these equations 

¥* to and m T can "be functions of time. This makes the integration 

’ L 

of these equations more complicated, "bat, when accomplished this 
would, through the relation "between e and <fi , yield the cyclic locus 
of- the shaft center, as well as the instantaneous resultant 
hydrodynamic force . 
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(i) Cyclic Squeeze Films in Journal Bearings 

In this section, the study is restricted to the particular 
case of dynamic loading, in which the journal does not rotate about 
its center hut moves under a variable unidirectional load. The load 
.vector is fixed in the space, and the shaft center follows a path 
determined by the nature of imposed load. 


Here , since 


03 = 03 ^ = o 


we have from eq.ns . (4*10) and (4,12) that 

5 

yecL R . . 2^ . 2 

w. M 6de 


t r 
2 


12 dt o h 3 f(H,£,h) 


and 


3 

pci, R 2 7r 2 Q 

I W cos 4 - • #- / 

0 12 dt o h 5 f(N,£,h) 


(4.13) 


(4.14) 


(4.15) 


Dividing, we get 


de 

d£ 


tan 4 


f x (e) 

e ^(iT 


where 

3 >) 

and 

f 2^ 


| 2 ^ sin 2 6 d8 
° h 5 f(H,A,h) 

j" 271 " cos^Q d6 
o h 5 f(N,A,h) 


(4.16) 


(4.17) 


(4.18) 
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Integrating (4*16), we get 


s in <j> = A 


exp [/ 


de 

e f 5 (e) 


1 


where A is the constant of integration and 


f j( £ ) 


^TeT 


(4 .19) 


( 4 . 20 ) 


Equation (4.19) describes the locus of the shaft centre 
under a unidirectional load, which depends upon the cost ant of 
integration A, determined hy the initial position of the journal. 

Equation (4*19) is independent of load which shows that the journal 
centre once located on a given path will always remain there regardless 
of the magnitude and nature of the applied load. The effect of 
load will only be upon the amplitude and speed on that fixed path. 

The substitution of <j> from eqn.(4.19) into either of the equations 
(4.14) and (4.15) would yield the relation between the load and 
eccentricity. 

(ii) Bearing Response to Sinusoidal Loading 

We now consider the case of shaft travelling along the 
vertical centre line(4>= o) , with an oscillatory load. 

W(t) = \ sin “pCi) (4.21) 

where W, is the amplitude of applied load and m is the frequency 

Jtr 


of oscillation 
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Equation ( 4 * 15 ) now "becomes 


77^ sinw^t 


ycL 5 E lc 2 7T 2 

o d£ f oos 9 de 


/ 


12 dt o h 5 f(lT,£,h) 


( 4 . 22 ) 


which after simplification becomes 

trc^y 


de 

dt 


12 \'L sinm t _ 

_b_ ■£_, 1 


( 4 . 23 ) 


uEcLq 2 


Using the non-dimensional scheme 


H = — = 1 +e cos 6 , L = 7 , T = a)' t = 27T nt 
c a p 


60 


720W 

p — i 

o lb 


( 4 . 24 ) 


and defining the Sommerfeld number, S, for an infinitely short 
bearing as 


q 2 _ dn Ax 2 
bl ur' - p ^c ' 

o 


( 4 . 25 ) 


eqn. (4.23) reduces to 


de 

dT 


sinT 


ttS (^) 2 f 9 (2ir) 


(4 . 26 ) 


where, D == 2R, and 


2w 2 


5 % / co s 6 d6 

f 2 ( 2 ir) -/ F(R,L,H) 


( 4 . 27 ) 


and 


E(H,L,H) 


12 


^ / i v2 _j__ ,, ra 

+ ^ “ 2LH C0th 2 


( 4 . 28 ) 



15 


Equations (4.26) and (4.27) are integrated numerically 

to yield the time variation of the eccentricity ratio over the 

region of T from o to 2it . The resulting response curves for 

various values of ITjL and S are shown in Pigs. 4*2 to 4.10. The 
L 2 

value of (-^2.) is taken to he 0.0 5. 

4 .4 RESULTS MP DISCUSSION 

One of the most important hearing characteristics is the 
response time, i.e. the time that will elapse for an oil film 
to he reduced to some minimum value, say a value which may have 
been established as the thinnest permissible film for the safe 
performance of a particular hearing. 

Pigs. 4.2 to 4.4 illustrate the journal centre locus 

for various combinations of the micropolar parameters U and L and 

the Sommerfeld number S, The value of (j~) 2 is kept fixed 
L 2 

((—.) = 0.05). Por smaller values of U and larger values of L 

the effects are not significant • . But for higher values of U 
and smaller values of L (where the micropolarity is significant) 
the effects are quite significant. This shows that the micropolarity 
causes the fluid to offer more resistance to journal motion, thereby 
allowing for smaller eccentricities (higher film thickness) for a 
constant load. 

The velocity of approach, de/dT, is plotted against time, 

T, in Pigs. 4. 5 to 4.7, over a half cycle. It can easily be seen 
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that the velocity of approach decreases for micropolar fluid for 
those values of T which are near to v/2, "but for those values of 
T which are close to 5 = o or T = t , the trend is almost same* 

The velocity of approach is highest for a Newtonian fluid "but 
for micropolar fluid this velocity decreases. Thus indicating 
an increased effective viscosity. 

The maximum eccentricity ratio, e , is a measure of minimum 
permissible film thickness, and thus its study is essential. The 
maximum eccentricity ratio is plotted against L in Fig. 4*8, against 
S in Figs. 4*9 and Fig. 4*10* In these figures it is seen that for 
larger values of N or smaller values of L, the maximum eccentricity 
ratio decreases. For the Newtonian case it is highest. It implies 
that the same amount of load is now being carried with a smaller 
eccentricity (or at a larger film thickness). This means that the 
system can be used to sustain much higher loads, or, conversely, 
the same load can be used for a longer period of time without any 
mal func t ioning . 



Eccentricity ratio 





4.3 Path of the journal centre o^er one cycle of sine 
loading for various values of coupling number 






Eccentricity ratio 








Fig 4. 7 Velocity of journal centre over one half cycle of sine 
loading for vari ous values of coupling number 




Maximum eccentricity ratio - 



L > 


Fig 4,8 Maximum eccentricity ratio vs. length ratio for various 
values of coupling number 







Fig, 4.9 Maximum eccentricity ratio Vs Sommerfetd number for 
various values of length ratio 
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CHAPTER 7 

LUBRICATION OP ROLLER BEARINGS III COMBINED 
ROLLING, SLIDING AND NORMAL MOTION 

5 .1 INTRODUCTION 

Thousands of papers on film lubrication have appeared 
since 1886 in the wake of Oshome Reynolds classical work on 
lubrication theory. In general, papers on film lubrication have 
extended the application of the fundamentals, as set forth by 
Reynolds, to a variety of bearing shapes with steady and transient 
films of variable density, viscosity and temperature. 

In recent years a lot of attention has been directed towards 
the study of rolling element bearings. Considerable understanding 
of the bearing characteristics has been obtained through the studies 
of elastohydro dynamic lubrication, and for all practical cases 

theoretical solutions have been published. 

* 

Most of earlier work in lubrication have dealt with steady 
lubricating films in which there existed either an effective 
tangential surface velocity or a maximum film velocity due to an 
externally pressurized source. Unfortunately there are few instances 
where the film are completely steady. Usually there are transient 
or periodic forces or displacements imposed on the films. Changes 
in surface velocity result in tangential accelerations which can 
be important during starting and stopping. How ever, relative normal 
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surface motions are of far greater importance. The film must maintain 
sufficiently precise spacing in the presence of disturbance. Various 
authors have examined theoretically the hearing behaviour under normal, 
sliding or rotating motion separately. Sasaki et.al [ 90 ] were 
the first who studied the dynamic behaviour of roller bearing. However 
the boundary condition applied at the trailing end violates the 
continuity of flux. Towson et . al. [ 91 ], later on, obtained a 
detailed theoretical solution of lightly loaded cylindrical bearings 
subjected to a combined rolling, sliding and normal motion applying 
a more realistic boundary condition at the film rupture point. A 
satisfactory agreement vra,s found between theory and experiment [ 92 ] , 

No attempt has so far been made to analyze this problem [ 91 ] 
from the microcontinuum view point . Thus it is felt that the present 
analysis may provide a better insight into the lubrication of such 
bearings when analyzed from microcontinuum point of view. Towson’ s 
et.al. results [91J are obtained as a special case of the present 
study. It will also be shown that this theory may yield better 
fit to the experimental curves obtained by Marldio and Towson [ 92 ] • 

5 « 2 THE PROBLEM ANT THEORETICAL -ANALYSIS 

The problem considered is that of the lubrication of lightly 
loaded cylinders in combined rolling, sliding and normal motion, 
with micropolar fluids. The fluid is assumed to be isoviscous. 

The geometry of the system and the coordinate system are shown in 
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Fig. 5.1. The two cylinders of radii Rj and R^, rolling with peripheral 
velocities and TJ^, respectively, and approaching with a normal 
velocity V , are separated hy a lubricant film of thickness h. The 
conventional representation of the system is a plane and a cylinder 
of radius R. Moreover, the circular cylinder can he approximated 
hy a parabolic cylinder with the film thickness h given hy [ 93 ] 

2 

, x 


(5.1) 


where 


1 _ 1_ 
R IL 


+ 


1_ 

=2 


where h Q is the minimum separation of the oil film between the 
surfaces . 


(5.2) 


Assuming the cylinders to he infinitely long in the axial 
direction, the governing equation for a two dimensional flow field, 
from eqn. (2.56), is 

~ {h 5 f(N,tyh) ~J} =p{ ~~ - 2 H + V} (5.5) 

where N and £ are given in (2.18) and f (ll,£,h) by eqn. (2.52). 
Integrating eqn. (5»3)» using the condition 

x 2 

" o» at X - x m or h - h m - h 0 + (5.4) 


the pressure gradient is obtained as 
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= 

dx 


V U 2 


lrf(N,ft f 'h) 


2 2 

(-jj-2) + 7(x- Xm ) } 


(5.5) 


Non-dimensionalizing, using 


X = 


(2Rh o ) 


w 


X = 
m 


m Jr h , . v 2 

"o' 1 ’ 


\ = 1 + x m » L 


, p = 


2 Pii; 


U (U 1 +U 2 ) f2W o 


1 ^2Eh 

4 ‘ WV 

dP = (x 2 4) + 2qQug 
^ h 3 f(n,l,h) 


(5.6) 

(5.7) 


where 


FO^L,!!) 


1 _J_ 

12 (lh) 2 ” 2lh 


coth 


HIE 

2 


(5.8) 


Equation (5.7) oan he written as 


dP _ 

ax ~ 


(X-X m ) (X+X m +2q) 

h 5 f(n,l,h) 


(5.9) 


Equation (5.9) shows that the pressure gradient falls to zero at 

X = X and X = -X ~2q. If the latter point is taken at the point 
m m 7 

of maxiraum pressure, then at the trailing end 


P «. 



at X = X 
m 


(5.10) 
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Condition (5.10) is the Swift -Stieber "boundary condition applicable 
when | cj < 1 [ 91 ] . 


The other boundary condition depends on the oil flow 
through the wedge. If the oil supply is enough and the weight of 
the oil is neglected, the pressure build up will begin at the start 
of the wedge, 'where X = _«, , This is the inlet boundary condition 

for the pressure, i.e. 


P = o at X s -00 (5.1l) 

Integrating eqn. (5.9) fear the pressure distribution, we 


have 


X (2q+X +X ) (XJC ) dX 

p(x) = / 2 


h 3 f(n,l,h) 


(5.12) 


The point X can be determined using P(X ) = 0, i.e. 
m m 7 

\ (2 q4X+X ) (X_X ) dX 

f m v nr x m' 

J - — =r » 0 (5.13) 

itf(it,l,h) 


Solving eqn. (5.13) we can find the values of X^ for different 
values of the parameters!",! and q. 


The load capacity per unit length is 


% a . 

W-/ p dx « - / 

«■ 00 


(5.14) 


or, the non-dimensional load capacity is 
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2%“ 


U (Ui + U 2 ) H/2Bh Q 


-- / 




x ax 3351 


(5.15) 


The shear stresses acting on the plane (y = o) and cylinder 
(y = h), in regions where oil is filled, are given "by eqn» (2.59) j 

hdE- 

T o,h ~ 2 dx + hg(ll,£ ,hj ' 


(5.16) 


where g(N,£,h) is defined hy eqn. (2.6l). In the cavitated region 
the mean shear stress is [94] 

T c = + ~ WfiTft 


(5.17) 


The surface tractions are 
*m 


dx 


h,2 ’ h ”I h 2 g (tr,a,h) 


( R .1 Pi) 


or, in non-dimensional form 



2h F, o 
0 


U (U 1 +U 2 )* / 2Hh 0 
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+ J 


r X m dX 

_L hg(h,l^)' 


H / 
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m 
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m 


dX 

H 2 G(h,L,H) _ 


(5.19) 


where 


G(h,L,H) 


1 


~ tanh 
hti 


MIB 

2 


1 ' Vh 


(5.20) 
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These surface tractions are put in the form 

F x = -Lb^+J C ? ] (5.21) 

f 2 = _ [B F „J C f 1 (5.22) 


(2l+X+X m ) (X^g d X 

^ -» 2H 2 F(1T,L,H) 

X 00 

f m <3X f dX 

F ~-L ^wi J H 2 0( . L H) 

m 

The coefficient of friction on the plane is 


( 5 . 23 ) 


( 5 . 24 ) 


h F 

(r) u f- (5 - 25) 

W 

The numerical solution consists of first solving the 
equation ( 5 . 13 ) to get the values of X^ for different parameters ; 

N, L, and q. The equation is solved using Regula falsi method 
within an appropriate accuracy where the infinite limit has been 
replaced by a finite but large quantity, (i.e. -» is replaced by -10) . 
A comparison is made between the solutions for the steady state case 
(q = o) with the non-steady case (q / 0 ) , For this purpose some 
ratios for different bearing characteristics are defined. 


(i) Location of the Film Rupture or Cavitation Point (X ) 
Let 

X 

x = m 


( 5 . 26 ) 
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where X^ represents the general point in which q can take any 

accepted values and (X m )^_ 0 is the point of cavitation for sliding 

or rolling case. For the same minimum film thickness, X is a 

rfl R 

function of micropolar parameters IT and L, and the velocity 
parameter q» 

(id) Location of Line of Maximum Pressure 

It has "been already pointed out that if X represents the 

point of cavitation then the expression for determining point of 

maximum pressure X , say, can he found as 

m p 

X = - 2q . - X (5.2?) 

m u m " 

P 

and thus the nature of this line would he same • Defining 




(X 


m q=o 
P 


(5.28) 


where q = 0 represents the steady state case, 
function of N,L and q. 


X 


m pR 


will he a 


(iii) The Maximum Pressure 


Let 



P 

m 



q-o 


(5.29) 
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where is the value of maximum pressure under pure rolling 

or s liding case for the same instantaneous value of the minimum 
film thickness , i. e. h Q » This is also seen to he a function 
of N, L and q. 


( iv) Instantaneous Load Carrying Capacity and Surface Tractions 


f E " S /® 4 =o 


and 


- 

^ = ^Vo 



°F 

«*)q-o 


(5.30) 


(5.31) 


(5.32) 


where the term ( )^_ 0 represents the steady state case, then it 

F and fL are dependent on 
R 1 R 2 

N, L and q. 


can he seen that the functions , 


The numerical solutions have been given for the point of 
cavitation, the instantaneous load carrying capacity and surface 
tractions. To have a better insight of the present analysis the 
pressure distribution, load carrying capacity and coefficient of 
friction are presented in non-dimensional forms for different values 
of N, L and q and compared with the Newtonian case (N = o or L -*■«>) • 
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5.3 BSBULSS AITS DISCISSION 


(i) 


Dimensionless Parameters s 


-^part from the micro polar parameters N and L, one more 
important parameter is encountered in the present study, i,e. the 
velocity ratio parameter, q, defined hy 


g. = 


Y 


V U 2 



(5.55) 


which determines the qualitative behaviour of a roller bearing when 
subjected to a combined rolling, sliding and normal motion. Positive 
values of q indicate a situation when the surfaces are moving apart 
normally and negative values signify an approaching case. Numerical 
solutions have been presented for ~o.9 < q < 0.9, since, it has been 
pointed out by Dowson E 92-1 , that the cavitation occurs for | qj <1. 


(iii) Bearing Characteristics s 


Non-dimensional pressure distribution curves are shown in 

Nig. 5.2 for different values of velocity ratio parameter q and the 

micropolar parameter N for a fixed L. Micropolar curves are shown by 
dotted lines. Nrom these curves it is seen that the approaching normal 
velocity (q negative) results in a stronger pressure generation while 
the separation (q positive) leads to a weaker pressure. The 
micropolarity ofthe fluid increases the pressure at corresponding 
points as compared to the Newtonian fluid (N = o.o). 
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Son-dimensional load capacity W is plotted against the 
length parameter L for different coupling numbers S and "velocity ratio 
parameter q in Tig. 5»3» If is clearly seen that as L increases 
the load capacity decreases and its value is seen to approach the 
Newtonian one for different values of q. For the fixed value of 
L the load capacity is higher for larger values of N and it is 
least at N = o.o (which corresponds to the Newtonian case [ 91 ] 
for the corresponding value of q) . A similar conclusion can he 
drawn regarding the parameter q. For negative values of q the 
load capacity is higher and for positive values of q it is smaller. 
When the normal motion disappears (q = o) , we get the case of pure 
rolling and sliding. 

The variation of the load parameter ratio Wjj with the modulus 
of the dimensionless velocity ratio q is shown in Fig. 5 . 4 . for 
different values of the micropolar parameter N. This figure depicts 
that tho effects of micropolarity aro higher for higher values of |q| • 
The parameter in case of normal approach, increases as q and 

N increase and decreases as these decrease (Fig. 5 .4) • The figure 
also shows that the load capacity is higher when the surfaces approach 
each other and it increases as the rate of approach increases or the 
micropolar parameter N increases. In case of separation the phenomenon 
is reversed. 

*"• «■* 

The frictional drag parameters F^ and F^ , 


as a function of 
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5T f L and q. have "been drawn in Fig-. 5*5* -An interesting feature 

of these curves is the graph for F„ './hen the surfaces are either 

“ ' R 2 

in normal appj?oech or normal separation. It is seen that the 

parameter 1 -^ first decreases and then increases when the surfaces 

are moving apart for all values of micropolar parameters. But when 

the surfaces are in normal approach the function F^ increases 

2 

slowly. The use of micropolar fluid is seen to reduce this factor. 

The qualitative behaviour of the parameter F^ is the same as that 
of load parameter. 

The coefficient of friction (R/h is plotted against 
the length parameter L for different values of the coupling number 
II in Fig. 5,6 (q. = -0.5), Fig. 5.7 (q. = 0.0) and Fig. 5.R (q_ = 0.25) . 

Afr r TJ 2 ) 

The value of Ja fprr r " is taken to he 0.5 in all these three 

1 2 

cases. It is seen that the coefficient of friction reduces to the value 
predicted by the classical lubrication theory for all three limiting 
cases (L -*■ o, L+», and H+o), It is seen that as the coupling 
number h increases the departure of the coefficient of friction 
from the classical theory becomes more conspicuous . The most 
interesting aspect of these curves is the fact that the coefficient 
of friction decreases rather sharply as L increases from 0 to 
approximately 10, The minimum value of is a function of L and 
IT which can be seen directly from the graphs. This observation is 
true for the case of normal approach and sliding or rolling. It can 
be also seen that as q. decreases from +0.25 (Fig, 5*8) to -0.5 (Fig, 5.6) 
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for fixed value of the coupling number N, the minimum values of 
coefficient of friction are obtained for smaller values of the 
length parameter L. Fig. 5*8 is the graph of coefficient of friction 

in case of separation of surfaces. A peculiar behaviour of coefficient 
of friction is seen for this case. As L increases the coefficient of 
friction first increases and attains a maximum value for a fixed 
coupling number and then decreasesrapidly and reaches a minimum 
values (a value which is always lesser than Newtonian one) and then 
increases again slowly. The behaviour is more significant for 
higher value of If, 

Fig. 5.9 is the graph for point of cavitation 3t m vs, the 
velocity parameter q., It is seen that as q. decreases the point of 
cavitation increases. The point of cavitation for pure rolling or 
sliding case is obtained at q. = 0 . 0 . In case of micropolar fluids 
the cavitation process is accelerated, i.e. it starts a little 


earlier 



4 cm ' 


q=0.50 


N = 0 6 


N=0.0 “\ 


' 7 / 


> 

" \ 1 





Fig. 5 2 Dimensionless pressure distribution 
of the dimensionless velocity ruco 
different coupling number N 













L 


Fig 5 7 Coefficient of friction (R/h 0 )^i f vs L for different 
values of N ( q =00) 
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Fig. 5.8 Coefficient of friction (R/h 0 ) pt vs. L for different 
values of N(q = + 0.25 ) 
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CHAPTER 71 

LUBRICATION THEORY FOR ROUGH SUREACiS AND 
ITS APPLICATION TO A JOURNAL BEARING 

6.1 INTRODUCTION 

The analysis of surface roughness effects on lubrication has 
gained considerable attention and has become a subject of intensive 
study during the last two decades. This is due to their importance 
in the practical field, its physical significance and also to the 
great theoretical interest inherent in these problems. 

Most theoretical studies on bearing lubrication appearing 
in the literature make use of the assumption that the bearing surfaces 
are perfectly smooth. It has long been recognized that owing to 
machining limitations, this assumption is rather unrealistic, 
particularly in bearings, where the film thickness is small. 

In recent years a major way of analyzing the problem has 
been the adoption of statistical methods. Two basic approaches 
have been developed namely the deterministic and stochastic approach. 

The starting point in each case has been the Reynolds equation. In 
the deterministic approach the surface roughness was accounted for 
by postulating a sine (cosine) wave or a series of sine (cosine) 
waves for the film thickness [ 95 J • More recently Tzeng and Saibel 
[96,97] have used stochastic concepts to study the surface roughness effects 
on slider bearing [ % ] and on short journal bearing [ 97 ] • The 
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analysis was concerned with, the problem of two dimensional slider 
bearing, with one dimensional roughness in the direction transverse 
to the sliding direction. There are, however, experimental evidences 
[ 98 ] which often exhibit a one-dimensional type of roughness 
running in the direction of sliding, generally referred to as the 
longitudinal roughness . Keeping in line with this evidence 
Christensen and Tonder [ 99-100 ] and Christensen [ 101] presented a 
stochastic model for hydrodynamic lubrication of rough surfaces . 
Christensen [ 101 J considered two types of one-dimensional roughness, 
the longitudinal roughness with its striation parallel to the sliding 
direction, and the transverse roughness with its striation perpendicular 
to the sliding direction, and obtained modified Reynolds equations by 
making two heuristic assumptions about the pressure gradient and flow 
rate. It was concluded that the longitudinal roughness results in a 
slight decrease in load carrying capacity and an increase in frictional 
force thereby causing a significant increase in coefficient of 
friction. The effect of transverse roughness is, however, to improve 
the bearing characteristics. Later Elrod [ 102 ] improved the 
Christensen’s modified form of Reynolds equation. 

' Scores of papers have appeared in the wake of the earlier 

mentioned works on various aspects of the lubrication of rough 
surfaces. All of the works appeared to date show that the surface 
roughness may have significant effect when the minimum film thickness 
is comparable to roughness mean height . The effect is much pronounced 
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for lower values of film thickness. Inherent in these works was a 
fundamental assumption that all the lubricants behaved as if they 
were Newtonian. 

Non-Newtonian behaviour is almost invariably observed in 
lubrication practice and they may be of important consequence when 
the surfaces are rough, Christensen’s models [ 101 ] do not account 
for non-Newtonian situations or where the micropolarity of the fluid 
is important. A question, therefore, naturally arises as to how 
the bearing performance would be modified in such situations. 

It is the aim of this Chapter to analyze the surface roughness 
effects from micr ocontinuum view point. For simplicity it is assumed 
that the moving surface is smooth and the stationary surface is rough. 
Three cases of roughness are considered and the theory is applied 
to the problem of journal bearing under half Sommerfeld boundary 
conditions . 

The present study may be of important consequence where the 
combined effects of surface roughness and additives are desired. 

6.2 MODIFIED REYNOLDS EQUATION FOR ROUGH SURFACES 

The Reynolds equation for micropolar fluids for smooth 
surfaces is of the form (eqn. 2.56) 

■h <v + ~iv ( v ■ tt" «•*> 



Ill 


where 


tih _ f ( k.j? „h) 9^ 

2. y 3x 


1 


z 


f(i T .ft.h) JP 
y az 


( 6 . 2 ) 

(6.3) 


l x an d ! z are obtained, from eqns . (2,5 3) and (2.54) by substituting , 

~ T Jj U 21 = ^12 = ^22 = °* ^ an ^ ^ are "by eqn. (2.18) 

and f(N,£,h) by eqn. (2,55), 


The geometry of the lubricant film can be considered to 
be constituted of two parts. The first part denotes the nomi na l (smooth) 
part of the geometry and can be considered to be a function of space 
and time co-ordinates. The second part of the film geometry, a 
randomly varying quantity, arises due to the surface roughness and 
is measured from the nominal level. Mathematically, 

h = h n (x,z,t) + h g (x,z,c) (6 A) 

£, a random variable, determines a definite roughness arrangement, 
from a large number of similar- but not identical - arrangements 
having the same statistical properties. 

Taking the expeoted value, equation (6.l) assumes the 

from 

~h tE( ^> + is {E W> - C6.5) 


where 
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£ ( ) - / ( ) £(h ) dh (6.6) 

o S 

— 00 

and f(h g ) is the probability density distribution of the random 
variable h g . Using (6.6) it is readily seen [ 101 ] , that 

^E(h) = 8h n 
at at 

where h^ is the smooth part of the film thickness* 

To evaluate the average of the flaxes appearing in equation 
(6.5), subject to a particular, specific model of the roughness, i.e. 
either along or perpendicular to the sliding direction, two assumptions, 
similar to those used by Christensen [ 101 j , are made . 

(i) The pressure gradient in the roughness direction is assumed 
to be a stochastic variable with zero (or negligible) 
variance 

(ii) The flux perpendicular to the roughness direction is assumed 
to be a stochastic variable with zero (or negligible) 
varianoe . 

A justification of these assumptions for micropolar fluid 
can be forwarded in a way similar to that of Christensen [ 101 ] . 

(i) longitudinal, One-Dimensional Roughness 

The roughness is assumed to be running parallel to the 
sliding direction in the form of long narrow ridges and valleys. 

Thus the film geometry assumes the form 


h = h n (x,z,t) + h s (z,5) 


( 6 . 8 ) 



113 


Taking the average of all terms in equation (6* 2) and adopting 
the first assumption the mean flux in the x-direotion is given "by 

= | E(h) - (6.9) 

This is so, "because, gp/^x is a variable with zero variance 
and f(lf,£,h) and Sp/s x can "be considered (approximately) to "be 
stochastically independent quantities. However, to get the mean 
flux in z-direction, we divide equation (6.3) "by f(H,£,h) . Taking 
the expected value of both sides of resulting equation and using 
the second assumption 

B(q ) » - ~ (6 .10) 

y 1 az 

E( ^f(N,a,h)' > 

Substituting the expressions for E(qp an( i E(l z ) ib- equation (6 .l) 
the modified Reynolds equation becomes 

-h- (E(f(lM,h)) 

- a- <n ^ * 2 -SsL) (6.11) 

fii) Transverse. One-Dimensional Roughness 

In this model, the roughness striation is assumed to be 
running perpendicular to the sliding direction in the form of long 
narrow ridges and furrows. The lilrn thickness therefore assumes 


the following form 
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h = h n (x,z,t) + h g (x,5) (6.12) 

In a similar way } adopted earlier, the flaxes, along and 
perpendicular to the sliding direction, are 


ECq*) 


(_ h 

u l 

2 tp/. 1 N 


1 3%(p) 

V 3X 


(E(f(lT,£,h)” 1 )) "- 1 


(6 . 13 ) 


E(l z ) - -i^|)-E(f(lI,£,h)) 

and the modified Reynolds equation is 


(6 .14) 



(iii) Uniform, Isotropic Roughness 

In this model, the roughness is assumed to he uniformly 
distributed over the hearing surface with no preferred position 
or direction on the surface. The film thickness assumes the 
form 

h = h^(x,z,t) + h s (x,z,g) (6.16) 

Adopting the same hypothesis as in [99 ] and mathematical reasoning 
given hy Tonder [ 103 ] , the fluxes and modified equations are 
given as follows . 
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S( V = f * y E (f(^»A»h) -^ L - 

E(q. z ) = -^ B(f(N,£,h) 


(6.17) 

( 6 . 18 ) 


is t E (fp.£.h)) yr 1 ) + -fr < fr 11 


JL. {u 9 g (h) + 2 5E (h) } 

2 3 x at 


(6.19) 


(iv) Distribution of Roughnc-ss Heights 

The roughness distribution function which is generally 
used to evaluate the several expected values is [ 101 ] 

f(h ) = -^ T (h 2 -h 3 ) 5 -b < h < b (6.20) 

3 32b 7 3 3 

= o elsewhere 

where h is the half total range of random film thickness variable 
and if a is the standard deviation, b = + 3° * The polynomial 

distribution function (6 .20) is an approximation to the Gaussian 
distribution. The reason for using such a polynomial distribution 
is that the Gaussian distribution always implies a finite probability 
of having asperity of very large sizes violating the conditions 
for hydrodynamic operation. 

6.3 BEARING CHARACTERISTICS 

Once the modified Reynolds type equation is obtained, for 
a particular case, the evaluation of the bearing characteristics 
can be carried out in a way similar to that of the smooth bearing. 
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(i) Oil Flow 


When the roughness is longitudinal the mean flux along 
and across the direction of sliding is given by equation ( 6 . 9 ) and 
(6.10), respectively. For transverse roughness these values are 
governed by ( 6 .13) and (6.1^). Equations (6 .17) and (6 .18) account 
for these characteristics for the isotropic case. Hie total flow 
over the edge of the bearing can be found from the above mentioned 
equations by integrating in a usual way. 

( id) Load Capacity 

The mean load carrying capacity is obtained by integrating 
the mean pressure over the bearing surface. 

E(f) = / E( p ) dA (6.21) 

A 

(iii) Frictional Drag 


The expression for shear stress on the moving surface, 
from eqn. ( 2 ^ 2 ), is 

T h 8P , uH 

T = 2 — + wTtir 


(6 . 22 ) 


where g(N,Ji,h) is defined by eqn. ( 2 . 64 ). 

For longitudinal roughness the mean shear stress is given 

■by 


( 6 . 23 ) 
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h \ 

fflff.&IhT ^ B (p) 

2_ 3 x 

'f (ii, a,h) A/ 


(6 .24) 


The mean frictional drag is 

EO) = / B(t) dA (6.25) 

A 

(iv) Coefficient of Friction 


The mean coefficient of friction is defined by 



( 6 . 26 ) 


6,4 MODIFIED RETHOLDS EQUATION APPLIED TO A ROUGH JOURNAL BEARING- 

The modified Reynolds equation derived above will now 
he applied to the lubrication of one -dimensional rough journal 
bearing under half Sommerfeld boundary conditions. 

The configuration of the system is shown in Fig. 6.l(<t) . 
Here the bearing is assumed to be rough and journal to be smooth. 
The two cases of roughness* i.e. longitudinal and transverse* are 


treated separately 
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(i) Longitudinal , One^imens ioual Roughnes s 

Since x = R$, equation (6*ll) assumes the form 

feMW)) SgLi-'f Mh 


de 


(6.27) 


The smooth part of the film thickness is 


h = c (1 + e c os 0) 


(6 * 28 ) 


and 

E(h) = h n (6.29) 

The boundary conditions for pressure are 

E(p) = o, at 0 =0 and ir (6'*3 0) 


The mean pressure distribution is given by 


0 ^ h n“ h n ^ d9 

E(p(e;) = pURj “E(f(fr }jl} h)) 


where 



ir h d0 

[ n 

J o E(f(N.a,h7) 

r" d0 

J 0 WWJM 


(6 «3l) 


(6.32) 


The mean load components per unit length, along and perpendicular 
to the line of centers , are obtained by integrating the mean 
pressure around the bearing surface from 6 = 0 to 0 « ir . Thus the 
mean load component normal to the line of centers per unit length 
is given by 
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IT 

EO^) = E(w) sin« = J 3(p(e)) sin.0 R d0 
2 ° 


(6.33) 


Substituting the value of E(p(@))from eqn. (6.3l) and integrating 
we get 


3(W*) 

2 


= | yI3R 2 / 


(h -h ) cob 0 d0 
n n 

o 


lE(f(lT,£,h)) 


(6.34) 


Similarly, the mean load component along the lino of centers 
is given 


E(W ) = E(W) cos <>= •“ pUR 2 / 


(h -h x . Q _ A 
tt n n ) sxn 0 de 


EU(N,i,h)')‘ 


The resultant mean load capacity is given by 

1 

E(W) - {(E(W 7f )) 2 + (E(W )) 2 } 2 


(6.35) 


(6 .36) 


The expression for mean shear stress assumes the form 
E (t) = + PUE( —j r-) 

J 2R de ^ ff.-M 1 


siW)' 


The mean frictional drag, e(p) , per unit length is 


given by 


- \ a" Ve® 1 E dfs +v ^ E( i5iy7hT ) E de (6 - 38) 


2 Tr 

J 

o 


(6.37) 


or, using eq.n. (6 .31) » 
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h (h -h ) d8 

E ( p ) = 4 ym{ i _ s{f(l,£,h)) ' 


+4 l E( 'iT!v;' 5 T ) 49 } (6 - 39) 


The coefficient of friction is obtained, according to 


equation (6 .26) . 


(ii) Transverse. One -Dimensional Roughness 


For this roughness the Reynolds equation is 




1 U TO t } (6 ,40) 

2 E <?edhf 


In a way similar to the case of longitudinal roughness, the mean 
value of pressure distribution, load carrying capacity, frictional 
drag and coefficient of friction are 


E(p(e)) = 


/ * E (f(N,£,h)^ ~ * d< 


n_^ ,h) 


(6.41) 


2^ / { E( f0^hT ) ” h n 1 E ^'f(lI,£, 


-) } cos 0 dg (6 . 42 ) 


:( w o ) -Jfa/ ^7^) ) -\ E( -W7v 


-)} sui 0 de 


(6 .43) 


L 0 El 'fp5 ,hj 




+ E( Rsytr } - 3 ~ •> 40 


E(77¥ 




(6.44) 
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where 


h 


n. 


/ E (' 


h 


-) de 


/ !( w’ d3 


(6.15) 


The coefficient of friction is obtained as usual* 


(iii) Isotropic Rotighness 

For one dimensional case the isotropic roughness will 
be equivalent to longitudinal roughness* 


6.5 hOH-DIMEHSIOiTAL FORMS 

To obtain the results numerically the following non-dimensional 


scheme is adopted, 
h 


h 


h 


h 


H 


s 


h = -r ~ 1 +£COS81 H o = r 2 * \ = r 1 * “s = v * 


B 


l ,£, s(p) , silk 2 , WJJ : aMg SWT - 


a UUR UUR“ pTJR 

(6*46) 

The non-dimensional average values of bearing characteristics 
for longitudinal and transverse roughness are given as follows. 


Longitudinal, one-dimensional roughness 


6 (H-H ) d9 

E(P(e)) - / s (l,l,h,b) 

0 -L 


( 6 . 47 ) 

— rzr rrTrrr . , f 

TT COS0d0 

'( 6 . 48 ) 

E(W 7r ) = E (¥) sin 4> = J 
2 


_TT 

(H-H ) sin 8 d 0 

( 6 . 49 ) 

e(W ) = E(w) cos d - / — 

0 n 

S 1 (lI,L,H,B) 
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, ir H(H-H ) de v 

2 J S 1 (N,L,H B) ‘ 2 / g G-S l » H » b ) d6 

0 -L 0 

(6.50) 

1 

E ('<?) = (s(Wp 2 + E(W o ) 2 ) 2 

(6.51) 

2 


_ m 

f(E(y f )) iM 

(6.52) 


where 


H <3.9 


E = 
o 


o 1 


s;(e,l,h,b)' 


TT 

/ 


cL0 




(6.55) 


S 1 (E,L,H,B) = 


+ 7 g ' (H + H s ) 2 ( B 2 -Ej) 5 coth f (H« ) 


G(E,L,H,B) = / 


B 


(B 2 -H 2 )d H 


(6 .54) 

(6.55) 


32B -B { ( H+H ^) _ Si tanh i^(H+E g )} 


Transverse, one -dimensional roughness 


E(P(9)) 

e 

= / {s 5 (e,l,h,b) - h o1 s 2 (h,l,h,b) > de 

(6 .56) 

i(w T 

TT 

2 

TT 

= J {S 5 (E,L,H,B) - H ol S 2 (e,L,H,B)> cos 0 de 

(6.57) 

E(W o ) 

7T 

= / {S 5 (E,L,H,B) - H ol S 2 (N,L,H,B)} sinede 

(6.58) 

W) 

if s 2 (u,l,e,b) 

2 { S 2 (N,L,H,B) “ H o1 S j( N s L » H » : b ) } d9 



2 IT S 2 (E,L,H,B) 

+ 2 J {3S 4 (E,L,H,B) +G(E,1,H,B) -3 -T^^l^bT 

> de (6.59) 
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where 


/ s (it,l,h,b) ae 

O J 


r 'ol f ir 


/ S 2 (3,L,H,B) d0 


( 6 . 60 ) 


S (N,L,H,B) - ^20 J 
32B' JB 


(H+H ) i “ 2 (B 2 JI ?) 3 a 


H 


[ (H+H g ) 5 + 12 - -^-(H+H s ) 2 coth ~(H+H s ) 1 

i=2 * 5 ' 4 ( 6 . 61 ) 

To determine the rate of change of hearing characteristics 

for a rough hearing as compared to the smooth theory the following 

hearing parameters are defined to study the micropolar effects. 

(iii) Mean Load Ratio Parameter 


The mean load ratio parameter is defined hy 

e(w J - (E(wyr B=o 


W R - 


( 6 . 62 ) 




B=o 


(iv) Mean Frictional Drag Ratio Parameter 


The mean frictional drag ratio parameter is defined hy 


S’. 


W) -Tmi 


B=o 


R 


( 6 . 63 ) 


(b(r)) 


B=o 


( v) Mean Coefficient of Friction Ratio Parameter 


The mean coefficient of friction ratio parameter is defined by 
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6.6 RESULTS AND DISCUSSION 

(i) Dimensionless Parameters 

Three dimensionless parameters are of important consequence 
in context with the present study, namely, N, L and B, defined 
by equations (2 ,18) and (6.46). 

The parameters IT and L have already been discussed in 
Chapter III. 

The parameter B( = “) is a characteristic of bearing 

c 

geometry and arises due to the surface roughness. It is obvious 
that the bearing characteristics would be prominent when b becomes 
comparable to c. 4 detailed discussion and limitation of this 
parameter are given in [ 99flOO,10l] • Here it suffices to say 

that o <_ B < 1. 

(ii) Limiting Cases 

The present study can be considered as a generalization 
of the works of Christensen [ 101 J and Prakash and Sinha [ 58 ] , 
because their results can be obtained from the present analysis by 
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considering certain limiting cases. 

Christensen' s results are obtained in the limit L » 
or IT -v o, as "both of these cases correspond to the Newtonian results . 
The results of Prakash and Siriha are obtained in the limit B -v o, 
since B o signifies a smooth surface. 

(iii) Bearing Characteristics 

The analysis of surface roughness effects from a microcontinuum 
view point, yields bearing characteristics which are functions, not 
only of the roughness parameter 3, but also of the micropolar 
parameters N and L, 

The load ratio parameter , friction ratio parameter 
and the coefficient of friotion ratio parameter C^ hare been plotted 
as a function of the roughness parameter B for various values of the 
coupling number N(e = 0 . 5 , L = 30) in Figs«6*2to 6 .4* It is seen 
that, akin to the Newtonian lubricants, the transverse roughness 
leads to a significant increase in load carrying capacity Pig. 6.2. 
and the frictional drag, Pig. 6*3 whereas the coefficient of friction 
decreases significantly, Pig. 6.4. The percentage rate of increase 
or decrease is higher for higher values of N. In case of longitudinal 
roughness there is a slight decrease in load carrying capacity, 

Pig. 6.2, and a slight increase in frictional drag, Pig. 6.3 
consequently the coefficient of friction increases for all values 
of N. Again it is seen that the rate of increase is further enhanced 
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"by large values of IT. 

Pigs .6 ,5 to 6.7 show the variation of the parameters 

F R an< ^ ^R versus 3 for various values of ll(e= 0.5 1 3=0.49) 

for longitudinal roughness. It is seen that the parameter W_ 

h. 

is lower for higher values of I and intermediate values ofL whereas 
is always higher. The reduction in W and increase in P leads 
to a significant increase in coefficient of friction and the rate 
of decrease or increase is always higher for micropolar fluid. 
However, the limits L -*■ o and L -*■ “ give rise to the Newtonian 
values for these parameters. 

f 

Pigs. 6.8 to 6.10 are the variations of the parameters 
F R and ^ for transverse roughness. and are positive 
showing that the load capacity and frictional drag are higher 
for rough hearing. However, negative values of show the 
reduction in coefficient of friction. The rate of increase of 
load capacity and frictional drag is higher for micropolar fluids. 
On the other hand the coefficient of friction parameter shows 
a typical trend. Prom Pig. 6.10 it can he concluded that for 
smaller values of L the rate of decrease of coefficient of friction 
(p f ) is smaller for micropolar fluids. However for higher values 
of L the decrease is higher and then comes to the Newtonian rate 
(corresponding to L ■> °°) . 
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Fig. 6-5 Load parameter Wr vs. length ratio parameter L for 
various values of coupling number N ( B=0'49 , 
longitudinal roughness ) 
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Fig. 6 -9 Frictional drag parameter F R vs length ratio parameter L 

for various values of coupling number N (B-0-49 transverse 
roughness ) 



•0 34 



833 1 52 24 36 48 L - > 60 72 84 96 
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of coupling number N (B = 0-49 transverse roughness) 
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CHAPTER 711 

SQUEEZE FIEfB BETWEEN ROUGH SURFACES 

7.1 I-TDROPHC TIOIT 

A number of experiments in which the evidence of rheological 
abnormalities, like enhancement of viscosity, has "been rather 
extensive and convincing, are of squeeze film type. For example 
Needs [ 22] conducted a series of experiments to detect the 
influence of boundary surface on the viscosity of lubricant films 
between two optically plane parallel circular disks approaching 
each other with no tangential velocity. He measured film thickness 
down to 0.000635 mm* For values less than 0,00127 mm an increasing 
discrepancy appeared between measured and theoretical (as predicted 
by the classical Newtonian theory) intervals of time. It actually 
took longer time for plates to approach than what was predicted 
by the theory, indicating that the effective viscosity of the 
lubricant in such thin films had somehow increased. Effective 
viscosity in the boundary film nearly five times the bulk value was 
observed. No conclusive explanation could be offered except to 
suggest that the proximity of the metal surfaces influenced the 
viscosity of the fluid, causing the fluid to become more rigid. 

Experiments by Fuks [28,29] also fall under the same 
category. Various fluids (like solution of pure fatty acids in 



139 


pure hydrocarbon solvent) were squeezed between two plates 
submerged in the fluid and measuring the gap between the plates as a 
function of time. It was observed that even after squeezing for 
a very long time a residual film remained between the plates. 

More recent experiments of this type are by Drauglis et.alL 31 ] 
and Askwith, Cameron and Gohar [ 36] ♦ In the experiment of Drauglis 
et.al. results similar to Fuks [ 28,29 ] were obtained, residual 
films were observed only in the presence of additives. However in 

Askwith et.al [ 36 ] squeeze film studies showed that the plates 

-4 

came to rest at a separation of about 2 x 10 ‘mm with all fluids, 
even pure Cetane. 

It has been suggested that asperities on the solid surface 
might account for the residual film . Also dirt and undissolved 
materials in the films have been invoked as the explanation for 
this residual film [ 30 ] • According to Drauglis et.al [31] 

"if the data of Fuks are all to be explained in this manner, then 
one must postulate dirt which changes its properties with temperature 
length of molecule of additive and concentration of additive but suoh 
dirt would be a most unusual material" • 

A proper theoretical explanation of the residual film reported 
by various workers is therefore clearly needed for a better understanding 
of thin film lubrication technology , especially when the height of the 
• roughness asperities is of the same order as ‘that of the thin films . 
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Motivated "by these considerations the squeeze film 
problems have been studied by various authors [ 55}68-71»75»76 J 
from aicrocontinuura view point . Inherent in all earlier theoretical 
analyses for squeeze films including those from the microcontinuum 
view point was a fundamental assumption of perfectly smooth surfaces. 

It has long been recognized that owing to machining limitations, 
this assumption is rather unrealistic, particularly in bearings 
where the bearing clearance is very small. Of interest, in this 
direction are the works of Prakash and Tonder [ IO 4 ] who studied 
squeeze films between rough circular plates and Prakash and 
Christensen [ 10 5 ] who studied squeeze films between rough 
rectangular plates. In analyzing the roughness effects use was 
made of the stochastic models for hydrodynamic lubrication of rough 
surfaces developed by Christensen and Tonder [ 99*100 ] and Christensen 
[ lOl] for Newtonian fluids. 

Hence in this Chapter, the lubrication theory for micro polar 
fluid for rough surfaces, advanced in Chapter 71, fs applied 
to the squeezing between parallel plates. Three standard geometries 
have been analyzed, i.e. two dimensional parallel plates (Fig.7*l(a) ) , 
circular plates and rectangular plates (Fig . 7 .1(b)) (three dimensional) . 
It is our strong belief that this theory could provide a proper 
theoretical explanation to the existence of residual films. 

The application of this theory may be found in the lubrication 
of rough surfaces where long chain length additives are used in 
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lubricants# An example to this context is the skidding of a rolling 
pneumatic tyre, even at low vehicle velocities on a Tret surface. 

The main feature contributing to this skidding is the existence of 
a squeeze film in the contact area between tyre and surface. This 
could be done by minimizing the- sinkage time. The behaviour of a 
tyre element in this zone may be idealized by the parallel plate 
squeeze film onto a rough surface and the fluid in between as a 
micropolar fluid on account of suspension of rigid particle additives 
(e.g. dirt particles) . A theoretical and experimental study of such 
a situation has been made by Moore [ 106 ] , where the asperity height 
of a randomly rough surface have all been assumed same. 

Further application of this study may be found in human 
joint lubrication where the mechanism is of squeeze film type. 

The lubricating synovial fluid in joints can be considered as a 
micropolar fluid and the cartilage surfaces rough ( a detailed 
justification is given in Chapter VIIl) . 


7.2 MATHEMATICAL ANALYSIS 

The generalized forms of Reynolds equation for rough surfaces 
from the micro continuum approach have been derived in Chapter VI. 

For squeezing, equations (6 .11), (6.15) and (6.19) reduce, respectively, 
to the following , 


Longitudinal, one-dimensional roughness, 
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{ E(f(N,£ ,h)) +-—{ (E(f(l?,£,hf 1 )-)" 1 -Jf^ } 


JEjh) 


3 

3 x 


3t 


Transverse, one -dimensional roughness, 


{(ECfdi,^)- 1 ))- 1 -^!} +~~ {E(f(:i,£,h)) —^1} 


a z 


32 


* v 


3E(h) 

at 


Isotropic roughness, 


{E(f(N,£,h)) +-|_ {E(f(h,£,h)) 


*p- 


3 z 

3E(h) 

3t 


where 


h = h n (x,z,t) + h a (x,z> 


(7.1) 


(7.2) 


(7.3) 


E ( )- / ( ) f(h s ) dh g 


f(h ) is the roughness distribution function of the random 
s 

variable h given by eqn. (6 .20). N and £ are defined in (2,18) 
s 

and f(N,£,h) by eqn. (2 .55) • 

Since the roughness distribution function is assumed 
to be symmetric, 


(7.4) 

(7.5) 


E(h) = h. 


n 


(7.6) 
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7.3 SQUEEZE FILMS BETWEEN ROUGH ITglSITBLY LONG PARALLEL HATES 

(i) Longitudinal , One-Dimensional Routines s 

For this type of roughness, the ridges and valleys are assumed 
to he running along x-direction. 

The film thickness assumes the form 

h = + \( z ) (7 *7) 

and the one-dimensional Reynolds eq.ua.tion, from eqn, (j .l) , is 

k - u ^ (7.8) 

The pressure boundary conditions are 

E(p) = o, at x = (7.9) 

MeI = 0 (e(p) maximum) , at x = o (7 .10) 


where is the half x-dimension of the plate. The mean pressure 
distribution is 


dh 


„/ \ 1 .. ~“n / 2 „ 2 n 1 


(7.11) 


The expected value of load carrying capacity per unit width is 

E(W) = 2 / E(p) dx (7.12) 

o 


or 


dh 


E(W) = 


1**3 
3 1 


n 


dt E(f 1 (S,£,h)) 


(7.13) 
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The non-dimensional load carrying capacity is 


E(W) = 


e (yi) 


> 3 dh 

/ 1 ru 

l ~3 dt j 


n 


P 1 (H f L,B) 


(7.14) 


where 


B 




■ / (B 2 -H^) 3 

32B _B 


( 1 +H ) 3 (1+H) 


12 


|;( 1 +H S ) 2 ooth } SH E 


(7.15) 


and 


h h . 

H - J , L - -S , B - f 
3 h n * h n 


( 7 . 16 ) 


The time of approach, t, from an initial (nominal) film thickness 
h^ to a final (nominal) film thickness h^» is given "by 




„ 3 h 

„ pi' n 
2 1 / 


dh 


n 


3 E(W)( 3(f (h,£ ,h) ) 

“i 

The non-dimensional time of approach is 

H dB 


T = 


1 ( 2 1 ■ \ 

h 2^ 3 E(w) ' 


f n n 

i 


( 7 . 17 ) 


(7.18) 


where 
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F 2 (lI, L ,H n ,B) = 


B _ (S+Hj 3 (H+HJ 

■XU_ J (b 2 ^ 2 ) 3 f 


32 B _B 


ns' ' n "s' 
*12 


n(H n +H s ) 


ooth ~ (H +H ) } dH 

9 > y \ a ' c 


2L 


(7.19) 


and. 


h h n. 

tj __ n ^ ^ s _ x ^ 1 

n ~ h ’ "“s h * a * ~ h 

n. n. n. 

xx x 


( 7 . 20 ) 


(ii) Transverse, One-Bimensional Roughness 

In this roughness, the ridges and valleys are assumed to 
Be running perpendicular to x-dir action. 

The corresponding Reynolds equation, from eqn » (7 «2 ) , 
is 




dh 


( 7 . 21 ) 


where 

h = h n (t) + h g (x) 

Proceeding in a way similar to longitudinal roughness 
the average value of squeeze film characteristics are given 
By 


9 * 

E(W ) 


-) 


Vi? \ 


¥T / e ( 7 Tw\" VT") 


t ' - 3 "WJJJ *4(1, £,h) J 

no. 


(7.22) 


(7.23) 


(7.24) 
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The dimensionless load capacity and squeeze time are 


E(Wj = f,(:t,l,b) 


H 

T = - / n F 4 (h,L,H n ,B) dH n 


(7.25) 

(7.26) 


where 


int - B (B 2 -H 2 ) 5 dH 

F,(N,L,B) = ^rr- J § s_ 

8B , 12(1+H ) 

{ (i+H y+ z-2- 


6N(l+H y 

j~-coth~(l+H g ) > 


(7.27) 


(B 2 -H 2 ) 5 dH 


P 4 (N,L,H n ,B) = 


B 

105 j “ 

8B? , 6N(H +h ) 2 ~ m 

{ (H n+ H s )hi2(B a ^ 3 )-- " .. a 


(7.28) 


and the other non-dimensional parameters for load capacity are 
defined according to eqn. (7.16) and for squeeze time as per 
eqn. (7 .20) . 

(iii) Isotropic Roughness 


The roughness is uniformly distributed over the surface 


of the plate. For two dimensional geometries the isotropic 
roughness will he equivalent to longitudinal roughness. 
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7.4 SQUEEZE FILM BEPWESU ROUGH CIRCULAR PLATES 
(i) Radial, One-Dimensional Roughness 

la this model, the roughness is assumed to have the form 
of long narrow ridges and valleys running in radial direction, and, 
therefore, film thickness assumes the form 

h = h n (t) + h g (e) 

The governing differential equation, from eqn. (7«l) in 
polar co-ordinates, is 


T^ } 


1 3 
“2 36 


3h 


(E (f (NjS-jh)” 1 ) ) -1 = 


For an axisymmetry case, eqn. (7 .30) reduces to 


dh 




The boundary conditions for average pressure are 
E(p) = o, at r = 

= o (E(p) maximum ) , at r = o 


where r^ is the radius of the circular plate, 


(7 .29) 


(7.30) 


(7.31) 


(7 .32) 


(7.33) 
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The mean pressure distribution is given by 


E(p) 


1L — Wr^-r 2 ) 

dt ^ 1 ' 


/ 


1 

E(f(sr,£,h)T 


(7.34) 


The average instantaneous load capacity and time of squeezing 
are 


J(W) = f 2ir r E(p) dr 
o 

y dh 


JL " ~n 4 1 

8 dt 1 ~E(f (N,£ ,h) ) 


t 


4 fr, 

tt p r J J 

" 8W h 


dh 

n 

E(f(N^,h)) 


The non-dimensional load carrying capacity and squeeze time 
are 


(7.35) 


(7.36) 


eIwT = 


E(¥) 


dh r!; 

/ 2LH E -i'i 

8 dt . y 

h 

n 


1 

^CN.LjB) 


t f* 

^ + 8 E(f T ^“2 


(7.37) 


(7.38) 


where P^(N,L,B) is defined by eqn. (7*15) and Fg^jLjH^B) by 
eqn. (7*19) and the non-dimensional parameters for load are 
according to (7.16) and for squeeze time as per eqn. (7 .20) . 
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(ii) Circumferential, One-Dimensional Roughness 

In this model, the ridges and valleys are assumed to he 
running along 6 direction. The film thickness is 

h = h (t) + h (r) (7.39) 

XI o 

The corresponding Reynolds equation, from eqn. (7.2), 
in polar co-ordinates, is 

ifj- i r(E(f(N,£,h)”3) -i ^1} 

+ -2-5T <E(f(lV,h)) -|§RR> (7.40) 

r 

For an axisymmetry case eqn. (7 .40 ) reduces to. 


{r(E((f( ■sj,,hf’)T 1 


= V 



(7.41) 


Proceeding in a similar way as in the case of radial roughness, 
the load capacity, squeeze time and their non-dimensional forms 
are given by 

_ ydh 

E ( w ) = - q r i E ^f ( n,a ,1^ 

Try 4 f „/ 1 

t " “ 8i W) r l 1 E( f(n/,h) 

n 
x 


) ah n 


(7 .42) 


(7.43) 
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and 


E(W) 


l&L 


( - 


* wdh n 4 
8 _dt 

n 


f 3 (u,l,b) 


H 


T = 


n 


TT ]i rr 


= - / F, (iT.L.H ,B) <3H 
' 4 * * n* n 


8h E(W) 

i 


where F 3 (N,L,B) is defined by eqn. (7 .27) and F, (h,L,H^,B) 
by eqn. ( 7 . 28 ) and the other non-dimensional parameters for 
load capacity are defined by eqn. (7*16) and for time height 
relation by eqn. (7.20) . 


(iii) Isotropic Roughness 

The modified form of Reynolds equation, from eqn. ( 7 .5) , in 
polar form, is 

-j^> *A,~- {B(f(jr I t,h))4^.> 


3r 


2ae 


ae 


a h 


=y 


n 


at 


For axisymmetry case the form of isotropic roughness will 
be identical to radial roughness. 

It is evident that the qualitative and quantitative 
effects of roughness will be same for both the geometries, i.e. 


(7*44) 


(7*45) 


(7*46) 
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infinitely long" parallel plates and circular plates (axisymmetry 
case). Hence the results are given only for circular plates. 

7.5 SQUEESIUG BETWBEET HOUGH RECTAUGUL.&R HATES 

The previous sections were concerned with the one-din ensional 
problems. Therefore the side leakage effects, the width/length 
ratio effects for unidirectional roughness were not detected. 

Moreover the distinction between the isotropic roughness and 
the one-dimensional roughness could not be made. Hence, in this 
section, the problem of squeezing between rectangular plates is 
analyzed in order to stui^r the effects of aspeot ratio (width/length 
ratio) on the lubrication of rough surfaces with micropolar fluids. 

The geometry of the system is shown in Pig.7*l(b). 

Since there is no relative tangential velocity between 
the surfaces and because of symmetry, the one-form of on e- dimensional 
roughness can be obtained from the other, simply by a rotation of 
co-ordinate axes. Hence, only the case of one -dimensional 
roughness, with roughness having forms of long narrow ridges 
and valleys running in the x-direction is studied. This will 
be termed as 'Unidirectional Roughness' • The case of isotropic roughness 
is also analyzed. 

(i) Unidirectional Roughness 

Since, in this model, the roughness is assumed to have the 
form of long, narrow ridges and valleys running in the x-direction, 
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the film thickness assumes the form 

h = h n (t) + h g (z) (7.47) 

The corresponding Reynolds type equation, from eqn. (7 .l) » is 




+ | rt »c*4,»o- J or l ^1) 


3 h 


= y 


n 

a t 


(7.48) 


Equation ( 7 .48) can be written as 


ifg&L + * 2 -ifg&L, 


ah- 


n 


ax 


R 


3 z 


E(f(N^ ,h))' 3 1 


(7.49) 


where 


6 R “ ~R 


6 p (R,il,h n ) = {E(f(R,£,h)) x S 


1 

2 


(7.50) 


Equation (7*49) is Poisson-type equation and can be solved by usual 
separation of variable techniques with the boundary conditions 


E(p) 

E(p) 


0 at x = 0 and x = 

\ 

o at z = + “ 2 **“ 


(7.51) 

(7.52) 


to give the mean pressure distribution 


E(p) = 


„ 2 dh 

4 y a ? n 

dt 


cosh 


— l 

,h)) n , 1>3i5 n 


n^ z 
a 2 6 R 


cosh 


nir b. 


- l)x 


srn 


nir x 


(7.55) 
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where is the x-dimension and hg 4 s 4he z -dimens ion of the 
rectangular plates. 


The mean instantaneous load carrying capacity is 
\/ 2 a 

2 

E(W) = / f E(p) dx dz 

-\A ° 


or 


16 V a? dh 

m - ~4 <tt) 5(f(K, t ,h)y 


/ “2 s H 

t — c tanh 


utt h 


h. 


n=l,5,5 n^ir 2a 2 6 R 2n 


The non-dimensional load carrying capacity is 
E(W) 


E(W7 - 


16 y A' 


n 


0 dh 

2 n 


dt 


(7.54) 


(7.55) 
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P 1 (N,L,B) 


co 6_ nirv. v 

I C Ep - tarJi “ .-i-) 
n=l, 5) 5 n tt 2a l 5 R 2n ^ 


(7.56) 


where 


6 g = fpg * A “ a 2 l3 2 » V A " lt~ 2 

P 1 (U,L,B) is defined hy eqn. (7.15) , P,(ir,I.,B) hy eqn. (7.27) 
and other non-dimensional parameters, i.e. B,L, H g are defined 


(7 .57) 
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according to eqn. ( 7 . 16 ) . 
The time of approach is 


48 U &r 
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The non-dimensional time of approach is 


T = 


t TT 4 eC’O h 2 
n i 


48 y A 


2a 2 6 R 


2n‘ 


4 ) dh n 


(7.58) 


+ v -2 / 


A H 


P 2 (i*,L,H n ,3) 


n 


niTV 

— s-tanh —“ 2 : 
n=l,3,5 tt xP 25 


I ( 


A 


R, 


2n 


,-)dH 

4 n 


where 


—2 

R 



(7.59) 


(7.60) 


The parameters H ,B,H and 1 are defined according to eqn. (7 .20) 
n s 

and R 2 (N,L,H n ,B) by eqn. (7 .19) and R 4 (R,L,H n ,B) hy eqn. (7.28). 
(ii) Isotropic Roughness 


The film thickness is 


h = h n (t) + R s (x,z,5) 


(7 .61) 
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and the Reynolds-type equation for this type of roughness is 
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}+ {E(f(u,£,h)) ^r^-} 
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(7.62) 


or 


3. E(p) + 3 . B(p) = 
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„ dh 
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The various bearing characteristics and their non-dimensional 
forms are given as 

Pressure distribution. 


4P a^ dh 


cosh 


utt z 

3 >r 


E(p) “ ^3 dt" 1 n-1,3,5^ Sin 
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(7.64) 


Load capacity, 


dh 


" K dt E(f(N,"£,h)) 


(7.65) 


where 
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K = I ( -f tanh ~ 
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( 7 . 66 ) 
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and is independent of the roughness parameter. 


Squeeze time, 


t = 


K 

kw 


n dh 

/ 


h E(f(N,£,h)) 
n. 

i 


(7.67) 


The load capacity and squeezing time can be put in non-dimensional 
form as 
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E(W) 


dh , 

K-^h 3 
dt n 




( 7 . 68 ) 


T = 
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= - / 
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dH 
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E(W)h‘ 


F 2 (H,L,H n ,B) 


( 7 . 69 ) 


n. 

1 


Thus it is evident that the qualitative effects of isotropic 
roughness on rectangular plate squeeze films are identical to that 
of radial one - dimensional roughness in circular plates (or 
longitudinal roughness in infinitely long parallel plates) « 


7.5, MCviERlGAL RESULTS 

The roughness effects on various bearing characteristics 
are studied through the following parameters 
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(i) Load Capacity Ratio Parameter 

The load ratio parameter, 1-^, is defined t; 


SO?)- (3(7/))-, 


=o 


R 


dlwJ) 


B=o 


(7 -70) 


where "(S(W) )-g_ 0 is the load capacity for smooth surfaces with the 
same nominal film thickness. 


(ii) Time Ratio Parameter 

The time ratio parameter, T^, is defined by 


1 


1 - Mb-o 

(T Vo 


(7.71) 


where (l)-g__ 0 represents the squeeze time for smooth surfaces 
to reach to the same nominal film thickness . While calculating 
this ratio the load is assumed to be same for rough surfaces 
and smooth surfaces as well. 


7.7 RESULTS MD DISCUSSION 
(i) Dimens ionless Parameters 

In addition to the well known parameters for smooth squeeze 
films , lubricated with Newtonian fluids , four new parameters 
N, L, B and are introduced in the present study. 

The micropolar parameters N and L have been discussed in 
Chapter III and the roughness parameter B in Chapter 71. Here it 
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suffices to say that the effects of surface roughness would he 
pronounced for large values of B and N and snail values of L. 

b 2 

The non-dimensional parameter (f— ) is termed as aspect 

Jx ^2 

ratio. In case of unidirectional roughness the role of this 
parameter might he important because the value of will determine 
whether the nature of the unidirectional roughness is longitudinal or 
transverse. 

(ii) Bearing Characteristics 

Bigs. 7.2 to 7*4 are the graphs for circular plates 
squeeze film and Pigs, 7*5 to 7*7 l° r rectangular plates squeeze 
film. 

Further, it should he noted that the cases of radial 
and circumferential roughness for circular plates are equivalent 
to the cases of longitudinal and transverse roughness, respectively, 
for the infinitely long parallel plates. 

Fig. 7*2 is the graph of vs. B where the nominal 
film height is kept constant , It is seen that when B i o,3 
the load capacity for rough surfaces does not differ significantly 
as compared to smooth ones. Even for micropolar fluid this 
assertion is true. Beyond B ► o.3 the roughness effects become 
prominent and as B + 1 the load capacity is affected significantly. 
It is seen that the radial roughness decreases the load capacity 
while, circumferential roughness increases it. The load capacity 
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for circumferential roughness may he times and for radial roughness 

O 

times the load for corresponding smooth surfaces, as S + 1, 

These limits are for Newtonian fluid and can he obtained directly 
from eqn.s (7.44) and (7.37) > respectively by taking the limits 
B -> 1 and IT -v o. For micropolar fluids the rate of increase or 
decrease in load capacity is higher and the deviations from Newtonian 
rate are more prominent at higher values of roughness parameter B and 
higher values of N* 

Fig. 7.3 and Fig. 7.4 are the graphs for T R vs. non-dimensional 

film height for B = o.l and B = 0.4, respectively. From Fig. 7.3 

it is clear that the effects of roughness are prominent only when 

the roughness parameter is of the same order as the film thickness, 

"B 

i.e. when jj approaches unity. It is also evident that the circumfer - 
n 

ential roughness increases the sihkage time while the radial roughness 
decreases it ani the effects of circumferential roughness are more 
prominent than the radial roughness. The reason being same as that 
given for load capacity. 

Fig. 7,4 shows a similar trend as shown by Fig. 7.3, except 
for the accentuation of roughness effects • 

Fig. 7.5 is the graph between f R and B. It is seen that 
the roughness effects are significant only for B > 0.4 and for 
values of not around 1,5* For values of around 1,5 the 
roughness effects are significant only for B very close to unity. 
Further, it is seen that the load capacity for rough plates is higher 
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than that for smooth plates for smaller values of (i.e. v , < l.p) • 
For values of from v = 1,5 to = 2.o, approximately, the load 
capacity decreases first, reaches a minimum, and then increases. 

For values of higher than 2 .0 the load capacity decreases for 
almost all values of 3, A mathematical reason for ITewtonian fluids 
is given in the following (for micropolar fluids a similar reasoning 
can be advanced, however, the mathematics involved will be complicated 
and hence only numerical results are given in this case to illustrate 
a similar phenomenon) . 

The unidirectional roughness pattern, with roughness ridges 
and valleys running in x-direction, tends to increase the x-component 
of the flow, ECq^.) , whereas the flow in z-direction, E(q^) , is 
reduced, mathematically, 






3 n ax 


(7.72) 


and 


B(q^) 


1 32b 7 

iai 35 


[ 3(5b^-b 2 )(b 2 -hJ) log 


h +b 
n 


h «b 


n 


+ 2bh n (l5h 2 -13b 2 )]” 1 ( 7#73 ) 

It is obvious from eqn . (7 .7 2 ) that the flux in x-direction 
increases and reaches to a maximum value as b + h n (or B + l), which 
is times the corresponding flux for smooth surfaces . The flux 
in z-direction decreases as-j“( = B) increases and as 1)4 h^(or B -v l) 
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3 

it becomes — times the corresponding flux for smooth surfaces. 

Further, it must he noted that for a constant area hearing 
the component of the flow along the smaller dimension would 
he larger because most of the fluid would have a tendency to escape 
along smaller dimension. Thus if v A < 1, < q z (f or smooth surfaces) . 

Thus when v . > v (v £ 1.5* v may he termed as critical 
aspect ratio ) hut not more than 2 ,o (approximately) , the x-dimension 
is smaller than the z-dimension and the increase in IE is more 
than enough to offset the decrease in E(<^) . This remains true 
for smaller values of B, For higher values of B, on the other 
hand, the increase in E(q_^) will not he able to balance the 
decrease in E(q ) . However when v . is sufficiently large, 
i.e. more than 2.o (approximately), the effect of decrease in 
E(l z ) will not he more than the effect of increase in £(<1^.), for 
a constant area hearing. 

Also, the roughness effects are more pronounced for 
micropolar fluids, i.e. the rate of increase (for smaller values 
of v^) and the rate of decrease (for higher values of v ) are 
higher when micropolarity of the fluid is taken into consideration. 

Fig. 7.6 is the graph, T^ vs. nominal film height H 
keeping the roughness parameter fixed (B = o.2). It is seen that 
for smaller values of V A the uni-directional roughness increases the 
sihkage time, i.e . the behaviour is similar to that of transverse 
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roughness. It is seen, that around v a = 1.5 the roughness effects 
are not prominent . For considerably high values of v^, the 
uni-directional roughness decreases the sinkage time for all values 
of B . 

Fig. 7.7 is the time - height relation for B = 0.4. 

The trend is similar to Fig. 7*6 except for the accentuation of 
the roughness effects due to the micropolarity. 

For isotropic roughness the load capacity and sinkage 
time are always lower than the corresponding smooth ones. It is 
seen that after >_ 15, the unidirectional roughness and isotropic 
roughness yield approximately the same results. The case of isotropic 
roughness from uni-directional roughness can be derived by letting 
v A — . 

It should, however, be remembered that the comparisons 
here have been made on the basis of nominal film height • Therefore, 
if h^ is the minimum film thickness for smooth surfaces, h — b 
would be mi nimum film thickness for rough surfaces and hence the 
minimum height permissible in order to avoid a surface - to - surface 
contact would be much less than that in the case of smooth ones. 
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CHAPTER VIII 

LUBRICATION OP HD2ISN JOII7T3 


8.1 IUTRODUC T ICP 

Two centuries of study has failed to unravel completely 
the mysteries of the mechanism of human joints, although from the 
engineers point of view they represent one of the closest points 
of approach "between the inventions of nan and the work of nature. 

The synovial joints provided "by the nature in the human 
body to carry out the trouble free motion of one bone past another, 
have long been identified as bearing systems. With the exception 
of the joints of the middle ear and that between the lower jaw and 
the skull, all synovial joints are potentially weight-bearing, 
although some exercise this function more often than others. These 
joints function as excellent bearings in the biological conditions 
and generally support considerable load, while providing, low friction 
service over a long span of one’s life time. There, exists no mechanical 
device which yields as low a coefficient of friction as that found 
in the synovial joints. Thus a clear visualization of the synovial 
joints is likely not only to provide a guide-lino to improve the 
existing mechanical bearings but also help in the replacement of a 
faulty human joint. 


The human joint, of course, is an excessively complicated 
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system from an engineering view point and it is not possible to 
build as a mechanical replica in which all the features of the real 
thing are represented, nor indeed is it possible to device a 
mathematical description that is completely comprehensive. It is 
more practical and probably more logical to commence investigation 
in the opposite manner, i.e. to start with the most simple squeeze 
film and to introduce complicating features one at a time, until 
either a situation of sufficient complexity is reached or the engineering 
and computational aspects inhibit further progress. 

This fundamental approach has led to the development of 
various theories of human joint lubrication. Serious study of the 
lubrication mechanism in human joints dates back to almost half a 
century ago. MacConaill [ 107 ] suggested as early as 1932, that 
the synovial joints function as the hydrodynamically lubricated 
bearings. Chamey [ 108 ] challenged this view in 1959* Sis 
experimental evidence supported a boundary lubrication action, 

Chamey observed a chemical affinity of hyaluronic acid molecules 
with the cartilage surface. McCutchen [ 109-112 1 introduced the 
porous nature of the cartilage into the lubrication concept, when 
he proposed his "weeping luhricat ion" ideas. Dint enf as s [ 113 ] 
put forward the concept of elastohydrodynamic lubrication for 
synovial joints, which was later supported by Tanner [ 114 1 • An y 
possibility of boundary lubrication was ruled out by Dintenfass, 
while Tanner indicated a slight possibility of boundary lubrication 
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at the asperity contacts . 

However, it was only in 1967 Dowson [ 115 ] realized 

that one mode of lubrication is insufficient to explain the joint 
mechanism. He concluded that "the major lubrication mechanism 
would seem to be some form of elastohydro dynamic action determined 
by sliding or squeeze film-action between porous surfaces with 
boundary lubrication providing the surface protection in cases of 
severe loading and little movement" . A little later in I968, 

Walker et.al [ 116 ], suggested a new concept of joint lubrication ,i, e . 
"Boosted lubrication". According to this concept of lubrication, 
squeeze-film action leads to a concentration of hyaluronic acid- 
protein complex in the lubricant as a result of diffusion of water 
and low molecular weight substances through the porous cartilage, 
and the restricated gap between the approaching cartilage surfaces. 

The increased concentration of hyaluronic acid molecule would 
give rise to an increase in viscosity of the synovial fluid in 
accordanco with the findings of ITegami [ 117 ] and the proposal 
is also consistent with the formation of gels on the cartilage 
surface (Maroudas [ 118 ] ) . Other work in this direction is by 
Higginson and Norman [ 119] who inferred that elastohydrodynamic 
lubrication plays a significant role, Unsworth ot ,al [ 120] and 
Wright and Dowson [ 121 ] have strongly suggested boundary lubrication. 

Though, Dintenfass [ 11 5} , Tanner [ 114 ] and Dowson [ 115] 
had made some numerical calculations for the fluid film thickness, 
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Fein [ 122 ] was probably the first to give a mathematical model for 
synovial joints . The joints system was considered to be elastically 
equivalent to sphere -on-a flat system. Mow and Ling [ 123 ] made 
use of Biot’s formulation of poroelastic theory for giving the 
displacement field equation and also the strain displacement 
relationship. The fluid flow in the porous matrix was taken to be 
governed by simple Darcy's law. Dowson et .al L 12 / 1-125 ] gave yet 
another simple formulation and accounted for viscosity variation 
with the hyaluronate concentration. The boosted squeeze time was 
found to be in good agreement with the experimental result. Ling 
[ 126 ] presented an. analytical model to account for the non-linearity 
of the cartilage. It was pointed out that the two mechanisms of 
•weeping', and 'boosted' are not exclusive. Higginson and Norman [l27 ] , 
Mans our and Mow [ 128,129 ] and Higginson [ 130 ] have given more 
rigorous mathematical models. They considered the elastic as well 
as the porous nature of the cartilage. 

It is however observed that the present situation regarding 
joint lubrication is not very clear. Earlier works based on the 
classical continuum theory, could not explain the observed effective 
increase in viscosity near the solid boundary. 

Little is known about the rheological nature and quantity 
of lubricant present in the human joints. Such knowledge is 
urgent ly required for the development of materials and designs and 
the better understanding of the tribological characteristics of 
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replacement "bearings for the human body. 

It would not be out of place to mention a little about 
the synovial fluid. Normal synovial fluid is generally clear or 
yellowish and viscous. It can be briefly described as a dialysate 
of blood plasma. It contains about one third of the protein 
concentration of the plasma. In addition the synovial fluid contains 
a very important polymer known as hyaluronic acid or sometimes called 
mucopolysaccharide. It is this hyaluronic acid molecule which gives 
synovial fluid its slippery and stringy behaviour. It also gives 
the characteristic non-Newtonian behaviour [ 131 ] . It was concluded 

by McCutchen [ 111 ] that the synovial fluid is a boundary lubricant . 

The synovial mucin are polymers of high molecular weight. Estimates 
of the molecular weight of the synovial mucin are as high as several 
million. Chemically, it is hyaluronic acid combined more or less 
firmly with less or more protein. 

Thus in studying the human joint lubrication problem, one 
must examine not only the continuum behaviour of the synovial 
fluid but also the microscopic structure. In obtaining the constitutive 
equations one must, therefore, resort to the microcontinuum approach 
as the continuum approach will be less definite about the detail 
structure of the material. 

This provides the motivation for the present work. The 
micropolar fluid model proposed by Eringen [ 44 I 




serves as a 
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satisfactory model for the description of the flow "behaviour of 
the polymer fluids and real fluid suspensions. It hr 3 already 
been pointed out that the micrcpolar fluid theoiy can he utilized 
to explain the rheological behaviour (increase in viscosity and 
boundary lubrication) of the lubricant. It also serves as a 
satisfactory model, when the fluid has long chain additive 
molecules . 

It is, therefore, noted that the synovial fluid may be 
thought of as a micropolar fluid because of the long chain hyaluronic acid 
molecules. The increase in the concentration of the hyaluronate 
in the synovial fluid, which is the important feature of the boosted 
lubrication, causing an increase in the effective viscosity, can 
also be accounted for in this theory, by way of increase in the 
micromotions and the couple stress. In this analysis, it is 
assumed that the human joint can be approximated by a system of a 
porous spherical bearing. Micropolar fluid theory is applied. 

Slip which almost invariably occurs near the porous boundary is 
considered by using the Beavers and Joseph [ 152 ] velocity slip 
criterion. The elasticity of the cart ilage ,i .e . the porous matrix, 
is ignored to avoid the additional complexity in the already sufficiently 
complex model. Thus the present analysis may be termed as a first 
approximation to the human joint lubrication theory, using the 
microcontinuum approach. The only work in this direction is by 
Tondon and Jaggi [ 75*74 1 • ^ho have used, rather too simple a 
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geometry for the human joint. The validity of various equations is 
debatable as vail be obvious from the present: analysis. The 
coeffic ient of friction is calculated wrongly, Moreover the full 
potential of the micropolar fluid theory has not been realised. 

8.2 MATOMT OP A SYHOVIAL JOINT 

The different forms of synovial joints in the Human body 
can not be readily reduced to a single model. For example, the 
hip joint is a close fitting ball and socket affair, to be 
approximated by an equivalent spherical bearing, while the ankle 
joint can be approximated by a cylinder-cylinder arrangement, 
resembling a close fitting journal bearing. Fig. 8.1, shows a 
joint illustrating the salient features encountered in synovial 
joints. In this schematic representation of a human joint, the 
load transmitting structural manbers are the bones. The ends of 
the bones which are usually globular in appearance, are covered with 
a layer of relatively soft and porous material called cartilage. 

The articular cavity contains the synovial fluid, which is 
responsible for the complex joint lubrication phenomenon. 

8.5 THEORETICAL FOMOLATIOH „ ASSUMPTION AFP EQUATIONS 

Geometrically the system of synovial joint, especially 
the hip joint may be represented by a sphere on a porous surface. 
Beth the sphere and porous surface are assumed to be rigid. The 
simplicity of the arrangement is shown in Fig, 8.2, The normal 




' -w 


Synovial membrane 
Articular cart lag e 


Fig 8 . 1 ( a ) Synovial joint 



<1 "'T , 

; t 



Fig 8.1(b) Hip joint 





180 


approach of the sphere is achieved by allowing it to fall under 
gravity on to the lower porous surface. The system is assumed 
to be lubricated by a highly non-lTewtonian 'Synovial Fluid 1 which 
has long chain hyaluronate molecules . liicropolar fluid theory is 
assumed to govern the behaviour of the synovial fluid. 


(i) Basic Equations 


The basic governing equations for a three dimensional 
flow field under the assumptions stated in Chapter II, from eqns. 
(2.32) to (2.35) ,are 


(v+ f*) 8 2u 
ay 


3v_ 


z +x ~ 


3 P 
3X 


= o 


( 8 . 1 ) 


(u+ ~x) — 2 -X' 
37 


3v^ 3p 

Ty " Tz 


= o 


( 8 . 2 ) 


3 2 v x 3 

__ _2Xv +X-22- 

3/ 1 


= o 


(8.3) 


3 2 v 


(8 .4) 


The boundary conditions are 


u - TJ n , w " V V 1 = V 3 = °» at 7 = ° 

/r 


(8.5) 


U = TJ, 
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3 u 

3 y 


y=h * W ~ U 22 
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a 3 ; 

3y 


y=h 


= o, = o, at y = h 


( 8 . 6 ) 
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where <j> is the porosity of the cartilage surface and A is the 

slip parameter. Condition (8*5) is the no-slip boundary condition 

for the velocity and microrotation velocity at the non-porous 

boundary and condition (8.6) is the slip boundary condition for 

to 

velocity according^ Beavers and Joseph [ 132 ] , however, the 

microrotation is assumed to vanish at the porous surface also. 

(ii) Velocity distribution 

Solving equations (8 .l) to (9.4) for velocity components 
u,w and microrotation velocity components v • with the boundary 
conditions (8.5) and (8 .6) j we obtain 

u m H f - B ny) - •? B ny - f 2 ( b 5 i sinB “y + B 41 cosh + B 21 

(8 .7) 

2 ^ jp ~ __2 

» - VS - ha y) - ~ h2 y - ~ (B 32 sM my + B 42 cosh my)+B 22 

(8 . 8 ) 

V 1 “ Tz - - B 32 cosh W ” B 42 sirih W (8.9) 

«*» T 3 

v 3 = " 2 rr( y Tx ~ ^i^ + B 3i 00sh ^ + B 4i sinl1 w (8.io) 

where 

" [ (hi-V + h i. { 241 i ^ * h - f 2 ^ f 1 ] 

{^(l-N 2 ) + h - ~ tanh —} 


(8 . 11 ) 
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7 -2 « p 

-> 0 . - - . ^ : — r — r — {B (cosh mh-l) + h — — } (8.12) 

2i xl p sinh mh. li^ 3 s. s 


B 0 • = B. .. + 


B,. = - ~±~ B, . 
3i 2P li 


B„. = 


4i 2u sinh mh li 


3p 

{ B.,.. (cosh mh-l) + h ) 


(8 .13) 

(8 .14) 


where 


1 

. X 2 jt 

i = 1,2, = x, ^2 - z, N = (~ + ‘x) > ® “ J » £ 


(— ) 2 
1 4y ; ’ 


(b.15) 


(iii) Generalized Reynolds Equation 

Integrating the equation of continuity 

(8 . 16 ) 

across the film, we get the generalized Reynolds equation as 


3 u 3v _ 3 _w 
3 x + 3 y 3 z 


3 P 


3p 


3 x { F il( lvI » £ » h » ^l) 3x } + 3 z { F ll^ n » £ » h » *1^ 3z } 


= 17 - U 11 (N,Sh,4 1 )}+-~ {G 12 (V,M 1 )> +V+7 


(8 .17) 


where V is the relative normal velocity of the surfaces and 1 is 
the velocity of the fluid towards the porous region at the cartilage 


surface and 
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f 5i { f n + ^ 1 ( 1 -- 2 )( 1-7^)} 
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G li- G li( E « s > h ) - 
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-L pi 


( 3 . 18 ) 


{ (l-N^h U u + | f 31 (lT li + U 2i ) } 

{^(1-K 2 ) + f 3] _} 


i = 1,2 

( 8 . 19 ) 

^ ^ r-.r „ , \ h^ , „ 2, IT £• h 2 I® 

f ll = = j2~ + 1 h ” 2 C0th 2£ 

(8 .20) 

f 21 ’ m . | + .“a 

sinh 

(8 .21) 

f 31 = “ k_2N£ tanh § 

(8 .22) 


Now the problem remains to find the Telocity of the 
fluid at the porous surface towards the porous region. For this 
purpose, let us recall that the hyaluronic acid i s a long chain 
polymer compound. The molecular weight of hyaluronic acid molecule 
is 5 94 10^ (approximately) and molecular length of the order 
5 x l(T 5 - 1CT 4 cm, lows on [115 ] • These molecules normally 
can not pass throu^a the porous matrix which is the basis of 
boosted lubrication. Thus it can be assumed that the flow 
through the cartilage matrix is Newtonian, and hence the velocity 
of the lubricant through the porous matrix shall be governed 
by the usual Darcy’s law [ 133 1 * 


Thus [ 133 ] , 
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V 
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6 

a 9P. [ 
y 3 y ! y=h 


( 8 . 23 ) 


Prom the requirement of continuity we have for the porous matrix 
$ ~ o 


e- 

V.q = - 


p = o 


(8 . 24 ) 


so that, since 



V 


2 

P = 


o 


( 8 . 25 ) 


The problem now reduces to get the solution of the Reynolds 

equation (8.17) for the pressure in the oil film, simultaneously 

with that of Laplace eqn. (8.25), for the porous matrix, with a 
3p 

common -r~ at the boundary. 

ay 

Assuming the thickness of porous matrix to be very small 
as compared to the radius of the sphere R) , the average 

pressure at any radial section of the porous natrix can be taken 
equal to the pressure in the lubricating film. It is further 
assumed that -Ir - is linear across the matrix and is zero at the 

ay 

outer surface of the porous bearing shell. 


This assumption leads to 
2 

— ^ ® constant = Kpsay) 

ay 


( 8 . 26 ) 


Hence from Laplace equation 
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Integrating (8.26) , we get 


DP 

37 


= K x y + K 2 


(8 . 28 ) 


where Kg Is constant of integration. 

„» 0 

With the boundary condition jf. « o at y - h + we get 


DP 

ay 


= K-, (y-h-H-j^) 


(8.29) 


Hence 


DP 


3yiy=h 


“ -hh 


Substituting the value of K^ from eq.n. (8.27), we get 


DP 

Tj 


. 3 2 p 3 2 
= H 1 ( -| + — 
y=h Dx 3 z 


(8 .50) 


(8 .51) 


Thus 


d> H 2 2 

P v 3x 2 Bz 


Substituting the value of V p in eqn. (8.17) and simplifying, we 
get 

]+ -ht t'uCH.i.Mii 


(8.52) 


9 H. 


l } _|£] = V + -~[{G 11 (H,a,h,9 1 )>3 


3z 




(8.55) 
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Equation (8.33) is the generalized Reynolds equation governing 
the pressure distribution in the cavity of the joint. 


8.4 H0EL1AL MOTION OF THE EDMAI7 JOUTO 

Attention is now restricted to the study of normal 
motion of the joints. Therefore, the effects arising cut of 
sliding are neglected, i.e. = o, i = 1,2. Thus eqn. 

(8 .33) takes the form 


3x 


[ { P 11 (H,^,h,6 1 ) + 


*a H l 


> 


+ -h [ + 


a 1 ^ 3 p -I 
U 3z J 


¥ 


( 8 . 34 ) 


The film thickness is 


h = c(l-e cos 9 ) 


and the velocity of squeezing is 

„ dh d £ 

¥ = — = -O TT" cos 8 
dt dt 


Using the spherical polar co-ordinates for the system under 
consideration, i.e. using 

x = R sin 9 cos t. 


y = R cos K 
z = R sin 6 sin C 
eqn. (s .34) assumes the form 


( 8 . 35 ) 


(8 . 36 ) 


( 8 . 37 ) 
(8 . 38 ) 
( 8 . 39 ) 
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(8 .40) 


It is reasonable to assume = o ( i . g . the pressure is constant 

3 £ 

along C direction). With this assumption, eon. (8 * 40 ) reduces 
to 

fr - [ sin e{P- n (N,£,h,cO + — 


r,2 . 39 

R sin 8 


'll' 5 * * V 1 
d£ 

= - c rr~ cos 0 
dt 


(8 41) 


To solve the Reynolds eqn. (8 «4l) the pressure boundary 
condition must be known. Since the pressure reaches a maximum 
at 8 = o, therefore the first condition becomes 


*E = 


d6 


o at 0 


(8 .42) 


However the other condition is debatable as the exact 
extension of lower cavity of bone is undetermined. If one assumes 

-X* 

the lower cavity to terminate at some points = + 8 f then the second 
condition would become 

p = o at 0 ~ + 9 (s *43) 

•5f 

But this again involves an unknown parameter 8 . However it is 
reasonable to assume that 9 = ~2 (approximately) . This assumption 
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may be justified because the pressure reaches a peak at 0 = o and 
falls very sharply on both sides, H e nce the influence of the 
boundary on pressure distribution is likely to be negligible 
after some distance from 0= o. 

Thus we have the second boundary condition as 
p = o at 0 = + 2 (8 ,44) 

Integrating eq.n, (8.4l) using the boundary conditions 
( 8 , 42 ) and (8*44) the equations governing the pressure gradient 
and pressure distribution, are, respectively 

IS - - H 2 o ff sine x 

„ ( ^(U 2 ) + f, x ) 

- hh * V + f 3i {r ii + ( i- _ 
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The expression for load capacity is 

2 //2 

W = 2tr R j p sin 9 cos d d6 
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(8 .47) 


or 
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Non-dim ensionalizing, using 
h 
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c = 1 -e cose, L = - 
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(8 49) 


where 4 and 4 are characterization of slip and porosity effects, 
P s 

respectively. Thus the non-dimensional load capacity is 
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where, 
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(8.51) 
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' l 2 L sinh 1TLH 
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2K , . 1TLH 

. ~ tanh — 
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The time of approach, t, from an initial film thickness to a 
final film thickness, in terms of corresponding eccentricity 
ratios, is given hy 
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The non-dimensional time of approach is 
2c 2 Wt 1 
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8*5 BESULTS AK1 3ESCU5SI0N 

(i) Dimensionless Parameters 

Apart from the usual dimensionless numbers encountered in 

the classical theory four new parameters $ , $ , IT and L, are of 

s p 

importance in the study of the present joint lubrication problem. 

The slip parameter <f> and the permeability (*> have since long 

s p 

been identified and discussed. It would suffice here to maition 

that a larger value of <J> or <j> signifies a greater slip or a 

s p 

greater permeability and a value of <f> g = o or $ = o signifies 

no slip or a non-porous body. It may be noted that for a non-porous 
body the slip is necessarily zero, however, a zero slip does not 
imply a non-porous case. 

The parameters H and L have already been discussed in 
Chapter III. In the present context, it suffices to say that, 1 
and D represent concentration and length parameter of the hyaluronic 
acid molecules, respectively, in synovial fluid. Higher is the 
value of N higher is the concentration of the molecules, and lower 
is the value of L higher is the chain length of molecules. Thus 
the joint lubrication characteristics would be prominent when N 
is higher and L is lower, 

(ii) Diseased Joints 

The properties of the synovial fluid undergo a marked variation 
with age, specially in case of a diseased or damaged joint. In a 
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rheumatoid arthritic patient, the synovial fluid may oven loose 

its non-Newtonian behaviour, Bloch and Dintenfass [ 134 ] • Bills is 

characterized by the cases of II = o or L •+ « , These cases actually 

correspond to the classical theory approach . Burch et . al [135 j 

indicated that the concentration of hyaluronic acid molecule in 

synovial fluid decreases for a diseased joint. For old joints, the 

cartilage becomes softer, less resilient and may even wear out in 

certain cases. These two cases can be characterized by a decreasing 

value of N and 9 respectively. Furthermore, it is seen that for a 
3 ? 

rheumatoid joint hyaluronic acid molecule has a lower molecular 
weight , A lower molecular weight corresponds to smaller length 
of the molecule. This, in this study, is depicted by a larger value 
of L. 

Thus it is evident that, the present study is capable, 

through the parameters <f> , <J> , L and N , to analyse the 

s p 

behaviour of a diseased joint also. 

(iii) Joint Characteristics 

The non-dimensional load capacity f is given by equation ( 8 .50) . 
Even a casual glance, tells that the slip parameter is multiplied 
by the factor (l-N 2 ), in the numerator as well as the denominator. 

The same is true for the expression of the dimensionless time of 
approach T, eqn. (8.55). This is an important result. This indicates 
that as N increases the effect of the slip will decrease, because 
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* 8 (u* ) will decrease-. Increase in IT, is a consequence of the 
Boosted lubrication. So as the concentration of the hyaluronic 
acid molecules increases, the effective viscosity becomes more and 
the effect of slip on the load and response time decreases. 

Curves for dimensionless load capacity *7 vs. the various 
parameters involved are shown in Figures 3 .3 to 8.9. Figures (8.3) 
to (8.5) show that the load decreases as the value of slip parameter 
increases. The qualitative trend Is similar to that obtained by 
the classical theory approach. Of importance is the fact that load 
decreases as (a) the porosity increases $ Fig. 8.3» Cb) the concentration 
decreases ; Fig. 8 .4 and (c) the chain length of hyaluronic acid 
molecule decreases ; Fig. 8.5* These results establish theoretically 
that a diseased or a damaged joint will carry a lesser load. 


Graphs for the load vs . $ ; Figs . 8 .6 - 8 »7» indicate 

P 

a similar trend that the load decreases as <f> increases for a 

P 

fioced value of other parameters. Fig. 8.7 5 which is a plot of 

W vs* 4> for different values of N, indicates that if d> 

P ’ P 

increases and h decreases, simultaneously, the decrease in load 
will be very rapid. 


Fig. 8^ is a graph of W vs. K for the porous (indicated 
by continuous lines) and the noruporous (dotted lines) cases for 
various values of L. The curve far L + ® gives the minimum load 
in both cases. This case corresponds either to the classical 



194 


theory approach or the reverse case of diseased joint ’."here the 
synovial fluid might have lost its non-iTewtonian characteristics. 
Obviously this curve is straight line, since it no longer is a 
function of II. Also, since the porosity of the cartilage has been 
established beyond doubt, the case of non-porous curves, in spite 
of their higher load capacity, are not of much significance . These 
curves are however given here only to illustrate the significant 
quantitative difference in the load capacity. 

Dimensionless load ¥ vs. L(Pig. 8.9) shows a very slow 
decrease in the load for low values of I, Indicating that for 
low concentration of- the hyaluronic acid molecules in the synovial 
fluid the load is very low, even when the chain length of the 
molecule may be high. But as the concentration further increases, 
the capacity becomes higher for long chain length (low value of l) 
of the molecule. However, as the chain length decreases, the load 
capacity decreases rapidly, irrespective of the value of the 
concentration. This happens in the diseased rheumatoid joint. 

The curves for the film thickness vs « the time of approach 
are shown in Pigs. 8.10-8.13. The study of these curves is very 
vital for the joint lubrication problem because this would indicate 
the time it would take the two cartilage surfaces to ccme into 
intimate contact, which actually will happen if the film thickness 
becomes zero, i.e* all the fluid is squeezed out from between the 


two surfaces 
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It is seen that the tine of approach is sufficiently large 

for low values of the slip parameter ; Pig. 8.10 , hut an increase 

in the slip decreases the time of approach ; thus preserving the 

qualitative trend of the classical theory approach. A similar 

result is visible from Pig. 8.11 ; that 1 decreases as the porosity 

increases. The time of approach for the slip case (continuous lines) 

and the no-slip case (dotted lines) is almost identical. The 

most important aspect of this figure is the fact that as the value 

of increases a finite time of approach is obtained. This becomes 

sufficiently small for high value of d j which indicates that all 

1 ? 

of the synovial fluid is either squeezed out or passes into the 
pores, causing an intimate contact of the cartilage surfaces in a 
very short period of time. This may be a cause of consistent severe 
pain in the diseased rheumatoid joints • 

Almost similar results are obvious from Pigs . 8 .12 - 8 .15 5 
wherein the results for the non-porous case are indicated by dotted 
lines and for that of the porous case by continuous lines. The time 
of approach is seen to decrease as either the chain length of the 
hyaluronic acid molecule decreases or as the concentration decreases 
(i.e. as L increases and H decreases) . Finite value of T are once 
again indicated for high value of L (low chain length) and low value 
of IT (low concentration) • 
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CHAPTER IX 


STEBI&Iff AMD COIICLUSIGSS 


I'he classical lubrication theory is based, upon a fundamental 
assumption that the fluid film behaves in a LTewtonian ray, i»c. 
on the classical continuum theory. However, the atomistic models 
of the material, especially of fluids that possess coarse structure 
and fibres (such as colloidal fluids, liquid crystals and fluid 
containing additives) have shown that the mass density can fluctuate 
violently with the size of the volume element, when it is below 
a certain critical limit but still in the macroscopic range. Thus 
the continuum assumption for mass density is no longer applicable. 

Another important class of problems where the classical 
theory is likely to give erroneous results is when the characteristic 
length scale is comparable to the average grain or molecule size contained 
in the medium. In such situations the molecular or the granular 
constituents of the medium are excited individually and the 
intrinsic motion of the mterial constituents must be taken into 
account. Typical of this class are the problem of flow through 
narrow passages, such as those found in human arteries and in 
lubrication theory. 

The classical lubrication theory also contains no mechanism 
to explain the presence of the residual film so commonly observed 
by various workers [ 22,28,29,31,36 ] . Thus there is a strong 
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motivation for extending the range of applicability of cent inarm 
mechanics, pa~*+icular]y in the s'hxdy of fluid flow of the type mentioned 
above, using tne field theory approach to describe the macroscopic 
manifestation of microscopic events, i.e, micromcti ons and defers at ions. 

Microscopic effects generated by the micromofc ions of particles 
in suspension of viscous fluids, or hy the molecules in the proximity 
of a solid surface, drastically change the character of flow between 
narrow passages. The origin of the rheological abnormalities is then 
to be investigated on the basis of the theory of micro-mechanical 
motions . 

One way of accounting for these micro-mechanical motions is 
through the so called theory of fluid microcontinua, advanced by 
Eringen [ 4-0 J , in which the continuous media are now regarded as 
sets of structured particles, which not only contain mass and 
velocity but also a substructure, i.e. each material volume element 
contains microvolume elements which can translate, rotate and 
deform independently of the motion of the macrovolume. Thus a 
mechanism is provided in this theory to treat materials which are 
capable of supporting local stress moments and body moments and in 
addition are influenced by the microelsnent spin inertia. 

The micro fluid theory is however too complicated and the underlying 
theoretical problem is not easily amenable to the solution of non-trivial 
problems in this field. This led Eringen to postulate a subclass of 
these microfluids, called micropolar fluids [ 44 3 , which would still 
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exhibit effects arising from particle micronotion. For this class 
of fluids, the deformation of the fluid microelement is ignored, 
nevertheless, microrotation effects are still present and surface 
and body couples are permitted. 

Physically the theory of micropolar fluids nay serve as a 
satisfactory model for description of the flow behaviour of polymeric 
fluids and fluid suspensions, or wherever the microrotation my 
play an important role. One is then naturally inclined to apply 
this concept of the microcontinuum to th-e analysis of a wide variety 
of fluid flow problems wherever a further refined model of fluid 
behaviour is sought . 

Motivated by these considerations the lubrication theory for 
micropolar fluids has been developed by many authors in the last decade. 
Various traditional lubrication problems have been reported. The 
results obtained are encouraging and illustrate the capabilities of 
this theory, in explaining the phenomenon of the enhancement of 
effective viscosity in thin films. Sven after a decade, this study 
seems to be in a preliminary state. Scores of lubrication problems 
have yet to be solved from this view point • For example, little work 
has appeared on dynamically loaded bearings and porous bearings from 
the microc csitinuum view point. Moreover, the problems of finite 
dimensional bearings and bearings with rough surfaces from this view- 
point axe yet to be analyzed. 
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Thus, ir. this thesis, these lubrication problems are studied from, 
the microc cut inuum point of vie'.?, in an attempt to enplain theoretically 
the rheological abnormalities so frequently observed in fluids with 
additiveo, especially when confined to narrow passages. 

The thesis begins with an introductory type Chapter, emphasizing 
the existence of the residual film in the proximity of a solid surface 
and laying the foundation for the application of the theory of micro polar 
fluids to ceitain hitherto untouched lubrication problems. Then in 
Chapter II, a generalized Reynolds equation applicable to finite lubrication 
problems, is derived, using the micropolar fluid theory. The equations 
obtained in Chapter II are applied in Chapter III, to a two-dimensional 
problem of squeeze film of a ball in a spherical seat and to some three 
dimensional non-cyclic squeeze films, assuming the characteristic 
coefficients to be constant, in an attempt to study the effects of rigid 
particle additives for the tnree dimensional micropolarity model. 

Increase in effective viscosity due to the micropolarity is established 
theoretically. It is also shown that the theoretical effects of the 
micro polarity on three dimensional lubrication are identical to the 
two dimensional problems, at least qualitatively. 

A very important result of this study is the dependence of 
bearing characteristics on a non-dimensional length parameter 

c Y 3/2 

L(= y , where c is the bearing clearance and £ = ( yyj-) ' , a fluid 
property having the dimensions of length) . This parameter links the 
fluid with the bearing geometry and can be thought of as a number which 
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is inversely proportional either to the length of the molecule of a 
non-ITewtonian lubricant or to the molecular sire of ~he additive. 

It is found that smaller the L(i«e. larger the H , or smaller the 
clearance) more pronounced are the micro polar effects. 

This dependence of hearing characteristics, especially effective 
viscosity, on the geometry of the system, may also give a possible 
explanation for discrepancies in actual viscosity measurements made 
using viscometers of different designs or sizes. 

It should also he noted that most of tie earlier workers have 
2 

used the relation o IS <1. However, in Chapter III it is shown that 
o 1 H 2 < -J • Thus the results of those workers are grossly exaggerated 
due to non-admissihle values of H, 

The next problem (Chapter 17) , deals with the dynamically loaded 
short journal hearings in which both wedge and squeeze films are 
active. The Reynolds equation for the general case of dynamically 
loaded infinitely short bearing is derived, where the lubricant 
is assumed to be micropolar* Detailed consideration is given to 
the dynamic behaviour of squeeze film In a short journal bearing 
under a sinusoidal load with no journal rotation. Various bearing 
characteristics are obtained assuming a full film to exist. The 
micropolarity of the fluid results in more resistance to journal 
motion, thereby allowing for smaller eccentricities for a constant 
load. The overall conclusion of this study is an increase in effective 
viscosity due to the micropolarity of the lubricant * 



2.12 


The next problem (Chapter V) , is that of the roller hearings 
in combined rolling, sliding and normal motion analysed under 
cavitation boundary.' conditions. Various hearing characteristics 
are obtained and represented graphically. The overall conclusion 
of the study, in conformity with the experimental results, is an 
increase in effective viscosity, resulting in an increase in load 
carrying capacity and a decrease in coefficient of friction, when 
analyzed from microcontinuum view point. 

It has been experimentally observed t 27 ^ that there is a 

sudden decrease in the coefficient of friction on addition of 

long-chain polar compounds . A graph of the values of coefficient 

of friction plotted against molecular weight (see ref. 156 , p.350) 

shows the greater the number of atoms and longer' the molecule the 

lower the coefficient of friction will be. Holding h constant 

° h 

(where h Q is the minimum film thickness) , and recalling that L = ~ 

(L is defined in this way in Chapter V) , L for the curves in Pigs . 5 . if 
to 5.g is seen therefore to depend only on the fluid, that is on 
the j material characteristic length. Considering the length of 
molecules of the boundary lubricant to be proportional to the material 
characteristic length, it is observed that the portion of the curve 
(N = 0 , 7 ) in Pigs. 5,£ to 5.g, to the right of the minimum satisfies 
the result that an increase in the length of molecule (corresponding v 
to a decrease in L) causes a decrease in coefficient of friction and 


v ice-versa 
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In the next problem (Chapter Tl) , micropolar fluid theory is 
applied to study the effects of surface roughness. Generalised 
forms of Leynolds equation are derived for three cases of roughness, 
namely, transverse roughness, longitudinal roughness and isotronic 
roughness, using the stochastic approach. For mathematical simplicity 
only the stationary part of the hearing - is assumed to he rough. This 
theory is subsequently applied to the problem of infinitely long 
journal bearing in order to study the interaction of micro polarity 
with surface asperities, using the half Sommerfeld boundary conditions. 
Various bearing characteristics are obtained numerically and results 
are compared with smooth bearing using nominal film thickness concept. 
The results obtained are in qualitative agreement with ilewtonian 
results for rough surfaces. However, it is shown that the effects 
of surface roughness are more pronounced, when analyzed from 
microcontinuum view point. 

A number of experiments in which the evidence of rheological 
abnormalities, like enhancement of viscosity, has been rather 
extensive and convincing, are of squeeze film type [ 22,28, 29f31>36] . 

It has been suggested that the asperities on the solid 
surface might account for the rheological abnormlities and residual 
film [ 30 ] . 

Motivated by these considerations, the next problem (Chapter VIl) 
deals with the squeezing between rough parallel plates * Three standard 
geometries, namely, infinitely long parallel plates, circular plates 
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and finite r-oiengular elates, are considered. Yard me r-ueezo film 
characteristics are- obtained for the three oases *f roughness . The 
results are c orrared with the correspond in." smooth plates. The 
rather important geometry, i.e. the finite rectan'ular places squeeze 
film is analyzed in details. Since this geometry can distinguish 
between the one-dimensional roughness and isotropic roughness, a 
comparision is made between these two cases. Further, for a constant 
area bearing, it is observed that the aspect ratio (width/length) 
plays a vital role for one-dimensional roughness, i.e. the bearing 
characteristics, for a given roughness, may increase or decrease 
depending upon this parameter. 

The last problem (Chapter Vlll) is of synovial joint 
lubrication studied from microcontinuum view point. The synovial 
joints provided by the nature in the human body to carry out the 
trouble free motion of one bone past another, have long been 
identified as bearing systems. The present situation regarding 
joint lubrication is, however, not very clear. Earlier work, 
based upon the classical continuum theory, could not give a 
satisfactory explanation for the increased effective viscosity 
in synovial joints. The microcontinuum approach represents a more 
realistic model for explaining the complex mechanism occuring in 
human joints. 

The behaviour of the synovial fluid, which has long chain 
hyaluronic acid molecules, is considered to be governed by the 
micropolar fluid theory. The porosity of the cartilage surface 
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and the slip at the porous boundary are tahen into account and a 
modified form of Reynolds equation is derived. The joints are 
approximated by a spherical bearing and various bearing characteristics 
are obtained in details . The salient feature of Lhis study, is an 
increase in effective viscosity. This theory also establishes 
theoretically, that for a diseased or a damaged joint the load 
capacity ana the response time decrease • This decrease may result 
in pain in case of diseased rheumatoid joint , 

The correlation of the results obtained in this thesis 
with various experimental studies may perhaps be believed coincidental 
and is definitely open to criticism since, in the light of these studies 
no experimental verification of the micropolar theory has been made. 
Nevertheless, at this early stage such criticisms are natural and in 
fact healthy for a crucial scrutiny of the theory and its ultimate 
acceptance. The results obtained are, however, interesting, and 
it is our strong belief that the micropolar fluid theory has the 
potentialities to explain the rheological anomalies in fluids 
containing additives through the parameters If and L when they flow 
through narrow passages. 

Although the various results based on this theory are in 
porfeot qualitative agreement with the experimentally observed results, 
an attempt for a quantitative fit would not be of any significance 
(since two undetermined parameters exist for curve fitting) unless some 
numerical values could be ascribed to the coefficient X and Y . 
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A significant step in this direction is by Arinar: et .cl [ 137 j 
who estahlished a relationship between tha ric_cpciar coefficients 
and hematocrit concentration in blood samples. Further, Vuglizrello 
and Sevilla [ 138 j , Goldsmith and fasen [ 139 ] , Sakhaveva 
et. al. [ 140] and ?uks [ 28,29 ] established a relationship between 
the effective viscosity and the molecular weight (particle size) • 

It is hoped that these studies would raise the curiosity 
of serious research workers to this rich field. The most useful 
contribution in this direction would be a repetition of the past 
experimental studies, in the light of these microcontinuum effects 
and thereby providing a much needed correlation with the theory. 
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